THE UNIVERSITY OF MICHIGAN 
GREAT LAKES RESEARCH DIVISION 



Technical Report 



NONDESTRUCTIVE ACTIVATION ANALYSIS! OF ENVIRONMENTAL SAMPLES 



Richard Dams 
John A. RolDldns 

John W. Winchester 
Project Director 



ORA Projects O89O3 and 125-^1 



supported by: 

U. S. ATOMIC ENERGY COMMISSION 
CHICAGO OPERATIONS OFFICE 
CONTRACT NO. AT(11-1 )-1705 
ARGONNE, ILLINOIS 

and 

INSTITUTE OF SCIENCE AND TECHNOLOGY 

THE UNIVERSITY OF MICHIGAN 

ANN ARBOR^ MICHIGAN 



administered through: 
OFFICE OF RESEARCH ADMINISTRATION ANN ARBOR 

May 1970 



TABLE OF CONTENTS 

Chapter I. Neutron Activation Analysis 

1. Principles 1 

2. Formation and Decay Rate of Radioactive Species 2 

3. Activation Analysis 3 

4. Sensitivity-Standard Deviations 4 

5. Counting-Gamma Spectrometry 5 

6. Nondestructive Activation Analysis . 8 

Chapter II. Procedure for Nondestructive Activation Analysis 
of Environmental Samples 

1. Apparatus Used 10 

2. Identification of Isotopes-Qualitative Analysis 11 

3. Quantitative Analysis of Germanivim Lithium Gamma-Ray 

Spectra 13 

4. Irradiation-Counting Scheme 15 

5. Sensitivity-Detection Limits 24 

6. Possible Interferences .... 26 

7. Neutron Flux Mapping 29 

Chapter III. Computerized Data Reduction 

1. Introduction 31 

2. Tape Analysis 

A. Tape Structure 34 

B. Tape Submission at U.M. Computing Center 34 

C. Inventory Program (INVEN) 35 

3. Data Analysis Programs 

A. Read and List Programs 39 

RDLST (main) 

RDTAPE 

MVL 

WRTMT 



B. Base Spectrum Calculation and Peak Integration . . 43 

BASE 
PKBS 
WRTBS 

C. Plotting 48 

D. Peak Search Routine 48 

PTST (main) 

PEAK 

SFZW 

FKZ 

CIJW 

E. Component Identification , . . . 51 

IDT ST (main) 

IDCOM 

GAMLIB (Data File) 

F. Non-linear Least Squares Program 55 

NLLS 
ARC 

4. Composite Program 59 

SPAN (Main Spectrum Analysis Program) 

III 

XM 

ECAL 

References 77 

List of Figures 78 



Appendix: A FORTRAN Listing of Programs 



INVEN 94 

RDLST 95 

RDTAPE 

MVL 

WRTMT 

BASE 

WRTBS 

PKBS 

PLOT 96 

PKTST 

PEAK 

CIJWZ 

SFWZ 

SZW 

FKZ 

IDTST 98 

IDCOM 
WRTID 
EF 

NLLS 99 

ARG 

SPAN 100 

XM 
IM 
ECAL 



I. Neutron Activation Analysis. 

Introduction 

The qualitative and quantitative determination of minute quantities by a 
technique known as activation analysis has become increasingly popular and is now 
generally recognized as a useful technique in various areas of scientific investi- 
gation and in an ever growing number of industrial control applications. 

In this technique a nonradioactive substance is made radioactive by particle 
bombardment, and its decay characteristics are then measured. The nature of these 
characteristics allow an identification and the magnitude is a measure of the 
concentration. The sensitivity is extremely high for a large number of elements. 
The specificity is usually excellent. A number of textbooks describing the prin- 
ciples, the possibilities and a number of applications are available (1-5 ). 
1. Principles i 

Activation analysis is the determination of the weight of a particular element 
in a sample by measuring the radioactivity induced by irradiation, usually with 
thermal neutrons, but in some cases with fast neutrons or even with photons or 
charged particles . 

The most commonly used and by far most intense neutron sources are provided 

by neutron chain reactors utilizing the fission reaction. Fast neutrons (>1 Mev) 

produced by the fission of uranium are moderated to epithermal ('^0.4 ev) and 

further to thermal (\0.4 ev) energies. These thermal neutrons assure continuation 

of the fission reaction but can also be used to produce artificial radioactivity 

induced when a target material is exposed to thermal energies a chance exists 

that another nuclide will be created by neutron capture, immediately followed 

by emission of electromagnetic radiation. 

The reaction is as follows 

A- A+1 
ZM (n,r) Z M 
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Examples are -i-iNa in,^) ^-.Na 

The nuclides formed, M, -j-iNa, 07^° ^^® often unstable and disintegrate, usually 
by emission of beta rays to an excited state of another isotope which deexcitates 
immediately by the emission of gamma rays. Beta rays are electrons or positrons. 
The beta rays from a given radioactive species are emitted with a continuous 
energy distribution extending from zero up to a maximum value. The gamma rays 
emitted by an excited state of a nuclide consist of electromagnetic radiations 
and have well defined energies E = hV (V is the frequency of the electromagnetic 
radiation). Gamma and beta rays can be detected with special counting equipment. 
More details about decay mode and interaction of radiation with matter can be 
found in the excellent books by Friedlander et al. "Nuclear and Radiochemistry" (6), 
and B. Harvey, " Introduction to Nuclear Physics and Chemistry" (7). 

2. Formation and decay rate of radioactive species. 

A target material A exposed to neutrons leads to the formation of nuclei B, 
which disintegrate to produce C as follows: 

A—^2k ^B — ^!^ ^C 

The rate of disintegration per atom of B is a constant, A, specific for nucleus 
B. Indeed within a fixed time (half-life = t 1/2) half of the radioactive nuclides 
B will disintegrate into C. ThusX= In 2/ t 1/2 (sec"-^) . The rate of formation 
of B depends on: (1) the number of atoms A present, namely N. ; (2) the neutron 

flux (j), the ntmiber of neutrons per second passing through an area of 1 square 

1 -2 -1 
centimeter (i.e. (j) cm -sec ); (3) the activation cross section of A, namely CT" , 

-24 2 
expressed in barns (1 barn =10 cm ) . 

The rate of accumulation of B is given by the rate of formation minus the 

rate of disintegration. 

Thus: ^=<r()N^-^N^ (1) 
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After integration assuming that before irradiation N ==0j one finds that at the 
end of irradiation t'^, the number of atoms of B present equals: 

N^i = HliiiA [1 _ exp (-\ti)] (2) 

The number of disintegrations per second, or the radioactivity at the end of 
irradiation t^ equals : 

Actgi = Ngi A = ry-^ N^ [1 - exp (- X t^)] = (T ^ N^ S (3) 
where S = 1 - exp ( -Xtj) is called the siaturation factor. The radioactivity 
present at time t^ after the end of irradiation will be smaller and is given as : 

Actgd = Actgi exp (-Xt^) (4) 

Thus knowing rj- , A or t 1/2, ^ and N the radioactivity for each irradiation 
time and cooling time can easily be calculated. Approximate values f or u'and 
t 1/2 are tabulated (6, 8). The neutron flux can be determined and is usually 

o 

approximately known for each irradiation site. Neutron fluxes vary between 10 

-2 -1 14 -2 -1 

n cm - sec up to several times 10 n cm - sec depending on the power of 

the nuclear reactor and on the irradiation site in the reactor. The number of 

atoms is given as 

_ 0.602 X 1024xgx9 ,^. 

A M.W. ^^ 

where g is the weight in grams 

M.W. is the molecular weight of the target material 
is the isotopic abundance 

3. Activation Analysis. 

Inversely one obtains the weight of the element from the activity measure- 
ment by substituting equation (5) in (3) 

Act^i X M.W. .,, 

g = B ^ (6) 

0.602 X 10^4 X ,rx ^ x S x 9 
This technique is called absolute activation analysis. The values for r are 
however not known accurately, a precise determination of (j) is time consuming and 
an absolute measurement of the radioactivity is difficult. Therefore, a comparator 
method is usually used in activation analysis. A standard containing a well known 



amoxmt of the element is irradiated simultaneously with and as near as possible 
to the unknown sample. A relative counting of both samples is performed and the 
unknown weight is calculated as follows : 

(g)v /(g) S = (Act)v /(Act)S (7) 

where (g)v and (g)S are the weights of unknown and standard (Act)v and (Act)S are 
countrates of unknown and standard. 

Many sources of errors can be eliminated by irradiating the sample with a standard 
of similar composition. When the composition of a substance is unknown, a prelim- 
inary irradiation has to be performed and can be used as a qualitative analysis in 
order to fabricate a suitable standard. Care should be taken that the sample and 
standard are approximately the same weight, shape and thickness. 
4. Sensitivity-Standard Deviations. 

Neutron activation analysis is inherently a very sensitive technique for 
determination of a large number of elements. In the analytical and radiochemical 
literature several mathematical expressions are used as definitions for the limit 
of detection, ranging from one to ten times the standard deviation of the back- 
ground and employing a widely varying terminology. Thus in order to draw valid 
conclusions with respect to detection capabilities from published data, it is 
necessary to examine carefully the definition and confidence level accepted by 
the author. 

As can be seen from equation (6) the sensitivity of the analysis of an element 
by neutron activation is determined by its cross section (a—), its isotopic 
abundance (0), the half-life of the Isotope produced (X), the Irradiation time 
(tj^), and the neutron flux ((|)) . Furthermore, the sensitivity is determined by 
the efficiency of the detector and its background signal. 

The following example shows how the sensitivity for an element can be calcu- 
lated from equation (6). 
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Assume the following realistic values: o^ =1 bam 

e =1 (100%) 

M.W. = 100 

ti = t 1/2 thus S = 0.5 

(p =10 n/cm -sec 
If the count rate necessary is 1 count per second with a detection efficiency of 

33 percent the sensitivity for such an element is given as: 

3x100 , ^_9 

S = 0.6 X 1024 X 2 X 1012 ^ 0.5 x 1 x 10-24 = 1° 8^^^ 

If gamma counting is performed the abundaii.ce of the gamma rays in the decay scheme 

must be included. 

Jenkins and Smales (9) calculated sensitivities assuming an irradiation at 

12 -2 -1 
10 n/cm -sec to saturation or up to one month (which ever is shorter) and 

counting 2 hours after irradiation to allow for chemical separation. The lower 

limits of precise measurements were taken to be 100 dis/min, for beta emitters 

and 100 photons /min for gamma emitters. Table 1 summarizes the results. 

As a matter of fact practical sensitivities may be much less favorable if 
a number of elements are determined simultaneously because irradiation cooling 
and counting time cannot be ideal for each isotope. An appropriate compromise 
must then be applied. 

Radioactive disintegration being a statistical phenomenon the standard devia- 
tion of a counting equals the square root of the number of counts. If the net 
number of counts N is obtained by subtracting the background B from the total 

number of counts T. The standard deviation is calculated as follows: 

'"2 2"* • — • 

^-n%-'t + ^^b %^^ + ^ 

5. Counting-Gamma Spectrometry. 

When a gamma ray strikes a scintillation counter (Nal(Tl)) it produces 

light flashes which release photoelectrons from a photosensitive electrode. 



-6- 



Element 


Estimated sensitivity, 
grams 


Eu, Au, Ho, 


10-12 


As, Ir, In, Kr, La, Mn, Pr, Re, Sra, 
Se, Tb, Tm, W, Yb 


10-11 


Sb, Ar, Br, Cs, Co, Cu, Er, Gd, Ga, 
Hf, Nd, Pd, P, Rb, Na, Ta, Th, U, Y 


10-10 


Ba, Cd, Ce, CI, Ge, I, Hg, Os, Pt, K, 
Ru, Ag, Te, Sn, Xe, Zn 


10-9 


Bi, Mo, Ni, Se, Si, S, Tl, 


10-8 


Ca, Cr, Fe, 0, Sr, Zr, 


10-7 


Pb 


10-6 



Table 1. Sensitivity of activation analysis as calculated by 
Jenkins and Smales (9). 



These are multiplied in a photomultiplier to produce a pulse which then can be 
further amplified and recorded. 

When a gamma ray strikes a semiconductor diode (Ge(Li)) a number of electrons 
are lifted from the valence band to the conduction band and swept through the 
material under the influence of an applied electric field so that an electric 
pulse is obtained. This pulse can also be amplified and recorded. 

Whenever the output from a detector-amplifier system is proportional to the 
energy dissipation in the detector the measurements of pulse heights may be very 
useful. Pulse height analysis is made versatile and rapid by the use of a very 
important tool in activation analysis, namely the multichannel analyzer. In the 
analogue to digital converter (A.D.C.) the pulses are sorted according to size 
and immediately recorded as counts in the appropriate channel of the memory unit. 
The resultant plot of the number of counts as a function of channel nimber is 
called a "gamma spectrum". The number of channels of an analyzer is determined 
by the capacity of the analyzer memory. Pulse height analysers with 50 up to 8192 
channels are in use. The time required for sorting and storing a pulse is appre- 
ciable and the analyzer cannot process a new pulse during this "dead time". 
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Provision is usually made to lengthen the counting time because of these dead 
time losses. A multichannel should also include the ability to display the gamma- 
spectrum on an oscilloscope and to drive an automatic printer, plotter or magnetic 
tape unit so that the content of each channel can be appropriately recorded. 

Gamma rays emitted by radionuclides being monoenergetic, it is possible to 
detect several isotopes in presence of each other if a response proportional to 
the original gamma-ray energy can be obtained from a detection unit. Thus if a 
scintillation crystal or a solid state detector is large enough to contain the 
entire path of the gamma ray, it can measure its total energy. This effect is 
called "photoelectric effect" and the resulting pulses corresponding to the total 
gamma energy form the so-called "photopeak". In addition to this there will be 
a number of "Compton" effects. This is an incomplete recovery of the gamma 
energy by the crystal and results in the appearance of a continuous distribution 
of smaller energies corresponding to the total energy minus the energy which escaped 
the crystal. This is called the "Compton continuum." High energy gamma rays 

striking a detector material can give rise to pair production. This involves the 

2 
production of a positron-electron pair (energy necessary is 2 mc = 1.02 Mev) . 

The positron is not stable and by the reverse process, positron-electron "annihi- 
lation", the same amount of energy is emitted as photons, namely two gamma rays 
with an energy of 0.511 Mev. If one of these or both escape the detector a peak 
is recorded with an energy equal to the gaiimia ray minus 0.511 or 1.022 Mev. They 
are called double or single escape peaks. 

The most important characteristic of the gamma spectrum is the full energy 
peak resulting from photoelectric absorption of the gamma ray. Thus the detection 
efficiency for a full energy peak and the ability to resolve it from another peak 
with a slightly different energy are of prime importance. The energy resolution 
is defined as the "full width at half maximum". Nal(Tl) detectors have high 
inherently full energy detection efficiencies and can be obtained in large sizes 



but unfortunately the energy resolution is rather poor (7.5% of the 0.662 Mev 

137 
gamma ray of Cs ) and is energy dependent. The Ge(Li) detectors on the other 

hand have lower intrinsic full energy detection efficiency and cannot yet be 

3 ■^ 

obtained in sizes larger than 100 cm . The detection efficiency of a 30 cm 

Ge(Li) detector may be 7 times lower at 0.2 Mev and 25 times lower at 1.3 Mev than 

the efficiency of a 3"x3" Nal(Tl) detector, but the energy resolution is much 

better, especially for high energy gamma rays. If coupled to good electronic 

apparatus an energy resolution of 2.5 Kev at 1.3 Mev is now easily obtainable. 

The resolution is only slightly energy dependent. Thus the resolution at 0.2 Mev 

is 8 times and at 1.3 Mev 25 times better than obtainable with Nal(Tl). 

More detailed information about gamma spectrometry can be found in the following 
works, "Gamma Spectrometry" by Heath, R.L. (10) and "Applied Gamma-ray Spectrometry" 
by C. Crouthamel, F. Adams, and R. Dams (11). 
6. Nondestructive Activation Analysis. 

The very high energy resolution obtainable in the Ge(Li) detectors allows 
the simultaneous specific measurement of the radiation emitted by several isotopes. 
Thus it becomes possible to distinguish and determine isotopes in complex mixtures 
without any chemical separation. The resultant instrumental or nondestructive 
activation analysis is attractive for different reasons. 

1) The sample remains essentially unchanged, and can be used for subsequent 
investigations after the decay of the induced activity. It also allows inter- 
laboratory comparison of results on the same specimen. 

2) Time is gained. This allows the analysis of large numbers of samples, 
which is very important in environmental research in order to draw valid conclusions. 

3) The instrumental approach allows the development of completely automated 
systems using computer handling of data. The analysis can be performed by an 
essentially unskilled operator. 



4) Isotopes with very short half-lives can be determined, whereas during chemical 
separations many short-lived isotopes decay essentially completely before the 
count is begun. 
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ll. Procedure for Non-Destructive Neutron Activation Analysis of Environmental 
Samples. 

Experimental 
1. Apparatus used. 

Nuclear Reactor: The irradiations are performed in the Ford Nuclear Reactor 
of the University of Michigan, operating at 2 Megawatts. The reactor is a swimming 
pool type, where the water acts as coolant, moderator and reflector. Short-lived 

isotopes are produced sufficiently by short irradiations (less than 30 minutes) 

12 2 
in the pneumatic tube facility (neutron- flux 2x10 n/cm - sec). The transit time 

of the sample is approximately 3 seconds and the sample holder is a special nylon 

tube (Rabbit). Longer lived isotopes are produced by longer irradiations (more 

than 1 hour) in a facility (H-5) in the immediate neighborhood of the core (neutron- 

13 2 
flux 1x10 n/cm -r sec). The samples are packed in a polyethylene bottle and are 

lowered in the pool by means of aluminum wire. Holes drilled in the side of the 

bottle provide sufficient cooling of the contents. 

Counting Equipment: The samples are counted by high resolution gamma-ray 

spectrometry. The gamma spectrometer consists of. 

Detector : Nuclear Diodes - True coaxial germanium- lithium drifted diode. It has 

a capacity of 27.5 pf. The detector is positioned in a vacuum chamber at 12 mm 

distance from a 0.5 mm aluminum window. Its performances are: a full width at 

fid 
half maximum (F.W.H.M.) of 2.5 Kev for the gamma transition of Co at 1332.3 Kev, 

and a detection efficiency for this photpeak of 4% relative to a 3"x3" Nal(Tl) 

scintillation detector. This results in a peak to Compton ratio of 18. 

Preamplifier : Nuclear Diodes-Model 101 is fixed to the detector chamber. 

Linear Amplifier : Tenelec Model 202 - an inverted unipolar pulse shape is used. 

It can be differentiated once or twice. 

Multichannel Analyzer : Nuclear Data 2200 Series - 4096 channel analyzer consisting 

of: Master Control; Analog to Digital Converter; Data reduce/integrate; Read in/out 

display; Memory Scaler; Memory Drivers; Memory array; Preset Counter; Magnetic 
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Tape Read; Magnetic Tape Control 

Oscilloscope ; Hewlett-Packard 

Printer ; Monroe-Data/Log - Model MC-13-80. It prints 20 channels per second. 

Magnetic Tape Recorder ; Kennedy 1400R - Incremental Magnetic Recorder 200 BPI. 

It is adapted for use with Nuclear Data 2200 Series. 

Counting Room ; The equipment is used in a temperature and humidity controlled room. 

Computer; IBM 360/67, with 2 CPU's sind other on-line accessories, at the 
University of Michigan Computing Center. 
2. Identification of Isotopes-Qualitative Analysis. 

High resolution gamma-ray spectrometry is very valuable for Identification 
of gamma-emitting isotopes. Indeed gamma-ray energies can be determined very 
precisely if proper care is taken to avoid instabilities and deviation from 
linearity of the analysing and amplifying equipment. Comparison of the gamma-ray 
energies with a compilation of energies of radionuclides formed by neutron irrad- 
iation (12 , 13) generally allows an immediate identification of the gamma emitting 
isotopes present without further investigation of other decay properties, or 
radiochemical behavior. Because the production reaction (s) of these isotopes in 
a reactor are also given in these compilations , the identity of the original target 
atom is usually very obvious. 

The spectrometer can be energy calibrated using a set of radionuclides with 
precisely known gamma- ray energies. Table 2 lists a number of calibration lines. 
For most applications a gain setting of 1 Kev/channel is very appropriate. Figure 1 
shows the deviation from this setting for the actually used gamma-ray spectrometer 
in an energy range from 55 up to 2800 Kev, a range which covers more than 95% the 
gamma rays emitted by reactor produced isotopes. By visual inspection of the 
oscilloscope display the location of sufficiently intense peaks can easily be 
estimated with an accuracy of better than 0.1 F.W.H.M. , which is here equivalent 
to approximately 0.25 Kev. 
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Isotope 


Energy 


(I 


Cev) 


241, 

Am 


59.568 


+ 


0.017 


131j 


80.166 


+ 


0.009 


^'\d 


97.43 


+ 


0.02 


^"cd 


103.18 


+ 


0.02 


^'hu 


112.97 






^^^Ce 


145.44 


+ 


0.05 


l^^Ce 


165.84 


+ 


0.03 


Lu 


208.36 






^'^Hg 


279.12 


+ 


0.05 


131^ 


364.47 


+ 


0.005 


198, 
Au 


411.776 


+ 


0.01 


Annihilation 


511.006 


+ 


0.02 


^O^Bi 


569.65 


+ 


0.10 


l^^Cs 


661.59 


+ 


0.07 


^\b 


765.83 


+ 


0.07 


^V 


834.84 


+ 


0.07 


88^ 


898.01 


+ 


0.07 


^O^Bi 


1063.82 


+ 


0.28 


6°Co 


1173.13 


+ 


0.04 


^°Co 


1332.39 


+ 


0.05 


Na 


1368.40 


+ 


0.04 


^°^Tl(ThC")(D.E.) 


1592.3 


± 


0.13 


^\a(D.E.) 


1731.6 


+ 


0.16 


88^ 


1836.1 


+ 


0.07 


^\a(S.E.) 


2242.6 


± 


0.14 


^°^Tl(ThC") 


2614.3 


+ 


0.09 


2'in.T 

Na 


2753.6 


+ 


0.12 



Table 2. Energy Calibration Lines 

From comparison with a calibration curve (Figure 1) an energy determination 
is immediately obtained with an accuracy of mostly better than 0.5 Kev. For more 
accurate energy determinations or for routine identification purposes a computer 
assisted qualitative analysis is performed. From a large number of calibration 
points a relation between gamma-ray energy and peak location is obtained and 
approximated by a nth degree polynomial. The degree of the polynomial and its 
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coefficients are calculated by a least squares fitting computer routine. Small 
differences in gain and zero threshold are corrected by recalibration with a 
high energy and a low energy gamma ray. The identification is then performed 
by comparing the experimental energies of the peaks of a complex gamma-ray spectrum 
with the gamma- ray energies of a library present in the memory of the computer. 
This library contains 2000 gamma-radiation energies of about 250 reactor produced 
isotopes (12). For the identification, also the half- life of the isotope is 
compared to the cooling time (time between end of irradiation and time of counting) 
of the sample. 

Chapter 3 about computer handling gives more detailed information about this 
qualitative analysis. 
3. Quantitative Analysis of Germanium Lithium Gamma- Ray Spectra. 

The quantitative determination of the count rate due to a certain isotope 

from the gamma-ray spectrum involves the measurement of the area included under 

its full energy peak. This determination must not necessarily be quantitative. 

It is indeed sufficient to measure only a fraction of this area, provided that 

the same fraction of the peak area in the standard sample is summed. This can 

easily be carried out as follows: (see Figure 2). If Ao is the count rate in 

the channel in which the highest number of counts is recorded, the area P under 

the peak is given as follows : 

n -n 

P = Ao + I^Ai + 'S.Ai 

1 -1 '" 

However the observed peak will usually be interfered by a contribution due to the 

Compton continuum of more highly energetic full energy peaks. The area Q due to 

the Compton continuum under the peak can be approximated as follows: 

^i^m -(2n+m+l) 
Q « Q' + 0" = -S-3 Ai + ^ Ai 
n+m - (n+m) 

The net peak area N is calculated as, 

N = P - Q 
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and bears (within the standard deviation of the counts obtained in the channels) 
a constant relationship to the total area contained in the peak. The method can 
be optimized in terms of experimental conditions such as the calibration and the 
resolution of the spectrometer and the values chosen for n and m. 

If the gain is set for 1 Kev/channel n is often taken equal to the F.W.H.M. 
of the peak, m is mostly chosen slightly higher. For the spectrometer used in 
this work the resolution varies from 2.2 to 3.0 Kev depending on the energy of 
the peak. During long counting times a small variation in gain, due to instability 
of the temperature in the counting room, is often experienced, resulting in an 
increased F.W.H.M. Therefore n can be chosen 2, 3 or 4 and for m a value of 4 
is mostly appropriate. The choice of n, depends on the stability of the gain, 
the complexity of the spectrum, the counting statistics and the energy of the 
peak. 

The precision of the determination of the peak area can be predicted if the 
peak is free of important interferences due to other peaks. The standard derivation 
is given as 





2 






2 






2 


Cf~ 


N 




-r 


P 


+ 


->» 


Q 


'J" 


N 


= 


J 


T 


'+ 


q' 



This method can however not always be applied: 

1) Sometimes another photopeak is present immediately above or below the 
photopeak to be measured, interfering with the calculation of Q' or Q". The 
Compton contribution is then calculated at another energy region which is free 
of interferences . 

2) If the photopeak itself is interfered by the presence of another photo- 
peak it is advisable to determine the abundance of the nuclide from another photo- 
peak. If this is not possible the interference may be reduced by decreasing n 

or the peak area may be calculated by doubling the area under the less interfered 
half of the photopeak. As a matter of fact the accuracy of such a determination 
is limited. 
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Because this nondestructive activation analysis procedure should be applied 
for large number of routine measurements, some sort of automatic data reduction 
is often necessary. Simple as these procedures for determination of count-rates 
under peaks may be, some sort of automation is sometimes desirable. Accuracy 
may be improved by the elimination of human errors in the transcription of data 
and in the operation of a desk calculator. However human judgment is to be main- 
tained by incorporating into the calculation procedures suitable criteria for 
checking data . 

A large number of methods and computer programs have been developed for the 
analysis of gamma-ray spectra obtained either with Nal(Tl) scintillation detectors 
or with Ge(Li) semiconductor detectors, but most of these programs involve thorough 
inspection and very complicated handling of data. Such a program is very expensive 
and thus less suitable for analysis of a very large number of quite similar samples, 
The next chapter describes the program for automatic data reduction applied in 
this work. 
4. Irradiation-Counting Scheme. 

High resolution gamma-ray spectrometry using germanium lithium drifted 

detectors allows the simultaneous measurement of a large number of isotopes. 

The ideal irradiation, cooling and counting times for the production and detection 

of an isotope are t, ^ t 1/2, t t *^t 1/2 and t ^ C^ 2 t 1/2. When a large 
^ irr ' cool' — count ^ 

number of elements have to be determined simultaneously by activation analysis, 
the irradiation-counting conditions can of course not be ideal for each isotope. 
Considering the elemental composition of typical environmental samples, the sen- 
sitivities for the determination of the elements present and the half-lives of the 
isotopes produced, a scheme, involving several irradiations and countings which 
allows the determination of a maximum number of elements was devised, 

a) Analysis of elements giving rise to short-lived isotopes: Preliminary short 
neutron irradiations of some environmental samples and several countings soon 
after the end of irradiation showed the presence of several short-lived isotopes 
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, 28, T 52,. 49„ 66^ 37„ 5L. 38„t 80„ 27^ 56„ 116m^ 
such as Al, V, Ca, Cu, S, Ti, CI, Br, Mg, Mn, In, and 

24 

Na. Table 3 shows the half lives and most prominent gamma- rays, for the analysis 

of the elements giving rise to these isotopes. The following irradiation-counting 
scheme was developed. Considering the half- lives of the isotopes produced it 
appeared appropriate to perform only a short irradiation and to start counting 
soon after the end of irradiation. A five-minute irradiation performed in the 

pneumatic tube from which a fast return of the sample is ensured. For such short 

12 2 
irradiation at a moderate neutron flux (2x10 n/cm -sec), the sample can easily 

be packed in a polethylene vial with a regular cap or with a snap top. Even 
■small amounts of liquid can be irradiated under these conditions. The vial 
containing the sample is placed in a special irradiation container (rabbit) and 
the experimenter himself can transfer his sample to the irradiation site through 
a pneumatic conveyor. Automatic or manual return after a predetermined time can 
be applied. Because the sample holder also becomes radioactive the sample is trans- 
ferred to another vial before counting. If the sample contains an appreciable 
amount of air it is advisable to ventilate it before counting because neutron 
irradiation of the argon present in the air (0.9%) gives rise to the production 
of the isotope Ar ^t 1/2 = 110 min) . For the same reason the irradiation container 
must be opened in a well ventilated hood. Three minutes after the end of irradiation 

a 400 second count is started. The bulk of the activity is then due to very 

.. , . ,- 28,, 52„ 49^ 5L. 66^ , 37„ ,^ ^, 

short-lived isotopes such as Al, V, Ca, Ti, Cu and b. Itoen these 

activities have died do-wn to a lower level a second count is started for 1000 

seconds at 15 minutes after the end of irradiation. This counting is favorable 

. , 38^, 27^^ 80^ 116m-^ 56„ , 24^, 

for the measurement of the isotopes CI, Mg, Br, In, Mn, and Na. 

Typical gamma spectra obtained by counting an environmental sample after 3 minutes 

and 15 minutes cooling are shown in Figure 4-5. The fixed irradiation-counting 

scheme is also illustrated in Figure 3. 
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Element 



Isotope 
Produced 



Half-life 



Most prominent Y rays 
Kev 



Al 

V 

Cu 

Mg 

Ca 

Na 

Ifo 

In 

CI 

Br 

I 

Ti 

S 



28 
52 
66 
27 
49 
24 
56, 



Al 

V 

Cu 

Mg 

Ca 

Na 



Mn 
116m 



In 



CI 



38 

Br 
128^ 



51 
37, 



Ti 



2 . 30 min 
3.77 min 
5 . 1 min 
9 . 5 min 
8 . 8 min 
15.0 hr 
2.58 hr 

54.0 min 
37.3 min 
18.0 min 
25.0 min 
5 .81 min 
5 . 1 min 



1778.9 

1434.4 

1039.0 

844.0; 1014.1 ; 

3083; 4071; 

1368.4 ; 2753.6 ; 

846.9 ; 1810.7 ; 2112.8; 
2522; 2657; 2962; 3371; 

138.4; 417.0 ; 818.8; 1097.1 ; 1293.4 

1642.0 ; 2166.8 ; 

617.0 ; 640.3; 665.7; 

442.7 ; 526.3; 743.5; 

320.0 ; 608.4; 928.5; 

3102.4 



704.3; 1256.7; 



Table 3. Nuclear properties of short-lived isotopes. The gamma-rays 
used in the determination are underlined. 



Standards: In order to perform a quantitative analysis it is necessary to 
relate the count-rate of each isotope to the amount of the target element present. 
Therefore it is appropriate to prepare a number of standard samples containing 
well known amounts of each element. These amounts have to be balanced carefully 
to obtain a reasonable activity of all the respective isotopes measured. 

Aqueous solutions were prepared from jmalytical grade salts of the elements 
to be determined. Table 4 shows the concentrations of the elements in these 
solutions. The solutions containing the metals consist of diluted nitric acid. 
The halogens on the other hand are present in diluted ammoniacal medium. From 
these solutions more diluted mixtures of the elements are prepared. Three mixtures 
respectively containing Al, V, Cu, Mg, Na, Ca, In, and Mn for solution 1; 01, Br 
and I for solution 2; and a separate solution 3 for titanium and manganese because 
the titanium solution contains an unknown amount of chloride and because the large 
amount of magnesium present in solution 1 interferes with the determination of 
manganese. Because these three solutions are very diluted with respect to several 
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elements it is advisable to prepare new solutions each 3 months. Amounts ranging 
from 50 to 200 X as shown in Table 4 of these three solutions are spotted on 
respectively 3 pieces of ashless filter paper (4x5 or 4x2.5 cm). Solutions 1 and 
3 are carefully dried under an infrared lamp while spotting. Due to the volalitity 
of the halogens, even in this ammoni a cal medium, this is not allowed for solution 2. 
These standard samples can be used several times provided that sufficient time is 
elapsed between two irradiations to allow for decay of the isotopes measured. 
Because sulphur is very insensitive by activation analysis a known amount of 
(NH,)„SO, can be weighed and irradiated. Amounts ranging from 50 to 200 mg are 
appropriate. 

The filter papers containing the standard solutions are then irradiated 
following the same scheme as the unknown samples. Because some of the elements 
spotted (Na, CI, Br, Ito.) are also present as impurities in the ashless filter 
paper, a blank paper has to be irradiated and the amounts found in this paper 
are summed with the known amounts to obtain the total standard. 

Flux monitor: Standard and sample cannot be counted simultaneously thus also 
the irradiation cannot be performed simultaneously. This assumes that the reactor 
flux does not vary significantly between the irradiation of standard and sample 
at the irradiation site. Frequent modifications in the composition of the reactor 
core makes this assumption however uncertain over a long period of time. Also 
the reproducibility of the irradiation position is unknown. Therefore a flux 
monitor is irradiated together with each sample and standard and used to correct 
for possible flux variations. A small foil of titanium (1 mg) is used as flux 

monitor and is counted at 13 minutes after the end of the irradiation for 20 seconds 

51 
The activity due to the isotope Ti (t 1/2 =5.8 min) produced by (n,y^ ) reaction 

is calculated and acts as a measure of the flux. This flux monitor is sealed in 

51 
a polyethylene jacket. Due to the short half life of the isotope Ti, the same 

monitor can be used almost indefinitely. Indeed, if the rate of analysis does not 
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Solution 


Element 


Concentration 




Volume 


Amoxmt 


Nimiber 




/ig/ml 


Spotted 
A 


Irradiated 


1 


Al 


10.0 


200 


2.00 


1 


V 


1.028 


200 


0.2056 


1 


Cu 


25.12^ 


200 


5.025 


1 


Mg 


500.5 


200 


100.1 


1 


Ca 


499.0 


200 


99.8 


1 


Na 


51.4 


200 


10.28 


1 


Mn 


1.214 


200 


0.2428 


1 


In 


0.076 


200 


0.0152 


2 


CI 


264.0 


50 


13.2 


2 


Br 


12,7-^ 


50 


0.636 


2 


I 


12. 8^ 


50 


0.64l5 


3 


r — ■ — -' 
Mn 


24.28 


100 


2.428 


3 


Ti 


415.0 


100 


41.5 


4 


S 


100 mg 
(NH4)2S04 


weighed 


24.2 mg S 



Table 4. Standard solutions oi: determination of elements 
giving rise to short-lived isotopes . 



exceed one analysis per 35 minutes more than 98% of the activity has decayed for 
the next measurement. 

Calculation of the Concentrations : The net-peak areas are calculated from 
the spectra obtained from standards unknown samples and flux monitor. The net 
activities induced in the standards and samples are divided by the activity in 
the simultaneously irradiated titanium flux monitor to obtain a relative activity 
per unit of integrated flux. The weights of the element (added and ashless 
filter paper blank) as present in the standards are then divided by the relative 
activities per unit of flux obtained for the respective isotopes. This yields 
a conversion factor for each element. This conversion factor is a constant and 
can be used for several months if the counting is performed in the same geometry 
conditions and if the performance and settings of the counting equipment do not 
vary appreciably. By multiplying the relative activities per unit of flux obtained 
from an unknown sample with this conversion factor, the total weight of the element 
(sample and blank) present in the sample is immediately obtained. By applying the 
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.fixed irradiation scheme no decay corrections are necessary because all the 

samples are counted at the same time after irradiation as the standards. 

The following formula shows this calculation. 

(m + b)i = (m + b)s • Is. • tL" K.fi 

As Fi Fi 

where: (m + b)^ = the weight of the element present in sample i (ff^^}^ "*") 

(m + b)„ = the weight of the element present in standard s (f ?° , ) 
^ blank 

Al = the peak area of the isotope in the sample i 

As = the peak area of the isotope in the standard s 

Fi = the peak area of the flux monitor irradiated with the sample i 

Fs •» the peak area of the flux monitor irradiated with the standard s 

K = the constant conversion factor for the element 

b) Analysis of Elements Giving Rise to Longer-Lived Isotopes: Because a 

long irradiation favors the production of longer lived isotopes, a sample was 

13 2 
irradiated in the pool of the reactor at a high neutron flux (1.5 x 10 n/cm -sec) 

for several hours. At different times after the end of the irradiation gamma 

spectra of the sample were taken. Up to 1 hour after irradiation the spectrum mostly 

shows the very short-lived isotopes discussed in the previous section. Between 

1 and 10 hours after irradiation the high activity due to Mn and Na greatly 

mask the presence of other isotopes. After 20 to 40 hours cooling time a number 

of photopeaks due to the following isotopes become visible: 

82^ 42^ 69m„ 64„ 122_, 140, 153^ 152m„ 24., , . 

Br, K, Zn, Cu, Sb, La, Sm, Eu, Na, and m some 

. 76, 72^ 187, 198, 
samples As, Ga, W, Au 

Isotopes with half-lives longer than 10 days only become visible in the 

gamma-ray spectrum after the complete decay of the isotope Br. Bromine is 

indeed usually abundant in environmental samples, it is very sensitive for acti- 
on 
vation analyses and the isotope Br has a very complicated ganma-ray spectrum. 

Twenty to thirty days after the end of irradiation the following isotopes can be 



-'^-^^-- .^^-^—---^—^—-^- 



21 



detected, ^^Cr. ^Sc, ^°Co, ^^e, ^^Zn. ^^Se. ^l%g, ^''s^, ^^^Ce, ^03 ^^, 



233 



58 



Pa. and in some samples Co, 



For the quantitative analysis of the elements K, Zn, Br, Cu, Sb, La, Sm, 



Eu, Ga, W, Au, Cr, Sc, Co, Fe, Se, Ag, Ce, Hg, Th, and Ni a number of environmental 



samples (up to 15) and a standard containing well-known amounts of these elements 

4 

irradiated simultaneously for 2 or 5 hours at a neutron flux of 1.5 x 10 n /cm -sec 




in the immediate neighborhood of the core of the reactor. The samples and the 



standard are heat sealed in polyethylene tubes. After irradiation the samples 



and standard are unpacked and transferred to the bottom of glass vials. After 



a cooling period of 20 to 30 hours the samples are counted for 2000 seconds for 



the measurement of the isotopes with half lives ranging from 8 to 50 hours. The 



isotopes with half-lives longer than 10 days are measured from a 4000 second count 



after 20 to 30 days cooling. Because the half- lives of the isotopes to be deter 



mined are sufficiently long as compared to the counting time a correction for the 



different decay times of the samples can easily be made. Table 5 shows the elements 



determined, the isotopes counted, their half lives and the energies of their gamma 



rays. Figures 7 and 8 show gamma spectra of both countings while the irradi 



procedure is illustrated schematically in Figure 6. 



Standards: Carefully weighed amounts of analytical grade salts of the 



elements are dissolved in dilute nitric acid. Well balanced amounts of these 



solutions are mixed to obtain reasonable activities of all the respective isotopes 



after irradiation. The solutions of gold contained some chlorine. Therefore it 



could not be mixed with the other elements (precipitation of silver and mercury) 



Also the addition of bromine and tungsten would cause evaporation (hydrob 



acid) or precipitation (tungstic acid H^WO , ) . Therefore these three elements 



are in a separate solution. The contents of these solutions and the concentrations 



of the elements are given in Table 6. Because the concentrations of several 



elements in these solutions are relatively low, new solutions should be prepared 
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Element 


Isotope 
Produced 


Half-life 


Most Prominent ^-rays (Kev) 


Na 


2^Na 


15.0 


hr 


1368.4; 2753.6; 




K 


^\ 


12.5 


hr 


1524.7; 




Zn 


69-Zn 


13.8 


hr 


438.7; 




Cu 
Br 


6^Cu 
^^Br 


12.8 
35.9 


hr 
hr 


f+ (511.0); 

554.3; 619.0; 698.3; 716.6; 827.8; 1043.9; 
1317.2; 1474.7; 




As 


7^As 


26.3 


hr 


559.2; 657.0; 767.5; 1228.8; 




Ga 


^^Ga 


14.3 


hr 


601.1; 630.1; 786.4; 810.5; 834.1; 1800.1; 
2201.4; 




La 


l^°La 


40.3 


hr 


328.6; 486.8; 815.4; 1595.4; 




Sm 


l^Sm 


47.1 


hr 


69.6; 103.2; 




Eu 


1^2-Eu 


9.35 


hr 


121.8; 344.2; 841.6; 963,5; 1315.0; 1388.9; 




Sb 


l^^Sb 


2.75 


da 


504.0; 692.5; 1256.6; 




W 


i«^w 


24.0 


hr 


72.3; 134.3; 479.3; 551.4; 618.1; 685.7; 772.9; 




Au 


l^«Au 


2.70 


da 


411.8; 




Sc 


^^Sc 


83.9 


da 


889.4; 1120.3; 




Cr 


5^Cr 


27.8 


da 


320.0; 




Co 


^«Co 


5.24 


yr 


1173.1; 1832.4; 




Fe 


^^Fe 


45.1 


da 


1098.6; 1291.5; 192.5; 143.0; 




Ni 


^«Co 


71.3 


da 


810.3; 863.5; f+ (511.0) 




Zn 
Se 


6^Zn 


245.0 
121.0 


da 
da 


1115.4; p+ (511.0); 

96.7; 121.8; 136.0; 264.6; 279.6; 400.7; 




Ag 


llOm^g 


260.0 


da 


657.8 


706.4; 763.9; 884.5; 937.2; 1384.0; 


! 


Sb 
Ce 

Hg 
Th 


^^^Sb 

l^^Ce 

^°^Hg 
233p, 


60.9 
32.5 
46.9 
27.0 


da 
da 
da 
da 


602.6 


645.7; 722.8; 1436.8; 1690.7; 2090.6 
311.8; 340.3; 




145.4 


279.1 


299.9 








._ 



Table 5. Nuclear properties of 
determinations are underlined. 



long-lived isotopes. The gamma rays used in the 
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detected. ^^Cr, ^^Sc, ^Oco, ^^Fe, ^^Zn, ^^Se, 1^°%, ^''s^, ^^he, ''\ and 
Pa^ and m some samples Co, 

For the quantitative analysis of the elements K, Zn, Br, Cu, Sb, La, Sm, 
Eu, Ga, W, Au, Cr, Sc, Co, Fe, Se, Ag, Ce, Hg, Th, and Ni a number of environmental 
samples (up to 15) and a standard containing well-known amounts of these elements 
are irradiated simultaneously for 2 or 5 hours at a neutron flux of 1.5 x 10 n/cm -^sec 
in the immediate neighborhood of the core of the reactor. The samples and the 
standard are heat sealed in polyethylene tubes. After irradiation the samples 
and standard are unpacked and transferred to the bottom of glass vials. After 
a cooling period of 20 to 30 hours the samples are counted for 2000 seconds for 
the measurement of the isotopes with half lives ranging from 8 to 50 hours. The 
isotopes with half-lives longer than 10 days are measured from a 4000 second count 
after 20 to 30 days cooling. Because the half- lives of the isotopes to be deter- 
mined are sufficiently long as compared to the counting time a correction for the 
different decay times of the samples can easily be made. Table 5 shows the elements 
determined, the isotopes counted, their half lives and the energies of their gamma 
rays. Figures 7 and 8 show gamma spectra of both countings while the irradiation 
procedure is illustrated schematically in Figure 6. 

Standards: Carefully weighed amounts of analytical grade salts of the 
elements are dissolved in dilute nitric acid. Well balanced amounts of these 
solutions are mixed to obtain reasonable activities of all the respective isotopes 
after irradiation. The solutions of gold contained some chlorine. Therefore it 
could not be mixed with the other elements (precipitation of silver and mercury). 
Also the addition of bromine and tungsten would cause evaporation (hydrobromic 
acid) or precipitation (tungstic acid H^WO, ) . Therefore these three elements 
are in a separate solution. The contents of these solutions and the concentrations 
of the elements are given in Table 6. Because the concentrations of several 
elements in these solutions are relatively low, new solutions should be prepared 
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Ag = the peak area of the Isotope in sample s 

'X = the decay constant of the element 

tj^ = the cooling time of the sample i 

tg = the cooling time of the standard s 
5. Sensitivity-Detection Limits 

In the irradiation counting system adopted a first detection limit is defined 
as the amount necessary to record at least 10 counts ((y = 35%) under the photo- 
peak or a number of counts equal to the background of the detector in that energy 
region. Practical detection limits are however mostly higher due to mutual inter- 
ference of the isotopes. Indeed the sensitivity to measure a photopeak of an 
isotope is decreased when the peak is sitting on the Compton continuum of higher 
energetic gamma rays "due to other isotopes . Thus a second detection limit may 
be defined which is defined as the amount of the element necessary to obtain at 
least a net photopeak activity equal to the standard deviation on its calculation. 
Thus if the net activity is given as N = P - Q, the detection limit is defined as 



v^- Q.' 

Column 5 on Tables 7 and 8 shows detection limits for typical air samples. 
As a matter of fact these values depend largely on the composition of the sample. 
For example a very high concentration of aluminum or chlorine may mask the presence 
of other short-lived isotopes . Large concentrations of sodium or bromine may 
interfere with the measurement of the isotopes with half-lives ranging from 8 to 
50 hours. This results in an important decrease of the sensitivity. 

The sensitivity of the analysis of an element in an environmental sample is 
further limited by the blanks. Indeed, abundant elements such as CI, Br, Na, K, 
Zn, Fe, Mn, Al, Cu, Sb, Cr, etc. are liable to be present in the filter material, 
holder, or solvent in which the sample is collected. 

In the present investigation an element is considered as being present if 
the net concentration exceeds twice the standard deviation of the blank value. 
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Element 


tirr. 
flux 


*- count 


decay 


Detection limit, 

ng, for typical 

air filter samples 


s 


5 mln 2xlOl2n/cm2-sec 


400 sec 


3 min 


25,000 


Ca 








II 


1,000 


Al 








II 


40 


V 








II 


1 


Cu 








II 


100 


Ti 








II 


200 


Br 




' 


1000 sec 


15 min 


20 


In 








" 


0.2 


Mh 








II 


3 i 


; Mg 








II 


3,000 '■ 


1 Na 








II 


200 


1 CI 








II 


500 


I 

! 








II 


• — 



Table 7. Sensitivity of 13 elements by short irradiation in adopted 
irradiation-counting scheme. 



Element 



K 

Zn 

Cu 

Br 

As 

Ga 

La 

Sm 

Eu 

Sb 

W 

Au 

Sc 

Cr 

Co 

Fe 

Nl 

Zn 

Se 

Ag 
Sb 
Ce 
Hg 
Th 



irr 
flux 



count 



5 hour 1.5x10-'- n/cm^ -sec 



2000 sec 



4000 



decay 



20-30 hours 



Detection limit, 

ng, for typical 

air filter samples 



sec 



20-30 days 



II 
II 
II 
II 
II 
11 



75 
200 
50 
25 
40 
10 

2 

0.05 

0.1 
30 

5 

1 

0.3 
20 

2 

1,500 

1,500 

100 

10 

10 

8 
20 
10 

3 



Table 8. 
s cheme . 



Sensitivity of 23 elements by long irradiations in adopted counting 
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This standard deviation is determined by the accuracy of the blank determination 
and by the homogeneity of the blank material. 
6. Possible Interferences 

Besides the thermal neutrons there are always some fast or fission neutrons 
present at each irradiation site of a reactor. These high energy neutrons induce 
threshold reactions a fter collisions with atoms. Such threshold reactions are 
(n,p), (n,c<), and (n, 2n) reactions. Some determinations based upon an (n,l(') 
reaction induced by thermal neutrons may be interfered by such threshold reactions. 

For example the determination of element pi by reaction (1) is interfered by 

A+1 A+4 A+2 
reactions (2), (3), and (4) on elements ^ M, „ ^j ^^'^ 7 ^^ 

^M (n,lf) f^ M (1) 

^11 M (n,p) f'^ M (2) 

Z+2 M (n,'X) 2 M (3) 

^+2 M (n, 2n) ^+^ M (4) 

To estimate the degree of interference it seemed appropriate to determine the 
ratio of the fast to thermal neutron flux. 

Determination of Fission to Thermal Flux Ratio: Determination of the fission 

flux in a reactor site requires the measurement of the production rate of an 

32 32 
isotope by a threshold reaction. Very often the reaction S(n,p) P is used to 

determine the value of an undisturbed fission flux. Thus from an amount of sulphur 

irradiated during a time t and counted at a time t^ after irradiation the fission 

flux is calculated as follows: 

= ^^32s exp (+Xt2) 



^r32s N 32g [1 - exp (-Xti)] 



/ 2 
where 0^ = the mean fission flux in n/cm -sec. 

32 
D^y = the absolute disintegration rate of P produced by (n,p) 

S reaction on sulphur 



g-^^ = the mean fission cross section for the reaction S(n,p) P 



N = the number of atoms S present 
J2g 
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X = the decay constant of the isotope P 

32 
This requires however an absolute measurement of the P disintegration rate. 

The interference due to a threshold reaction does not depend on the absolute 

value of the fission flux but only on its ratio to the thermal flux. Thus it is 

32 
sufficient to compare the count rate of P produced by the threshold reaction to 

32 31„ 32 

the count rate of P produced by the thermal reaction 'TCnj'Jf) P. This requires 

32 
now only a relative measurement of P. Thus from a simultaneously irradiated and 

counted amount of phosphorus the thermal neutron flux is given as: 

d = D31p exp (+Xt2) 

th 



(7- , ^22 ^-^ ~ ^^P (Xt^)] 



where 033^ = the absolute disintegration rate of P produced by (n,?r) reaction on 
phosphorus 

31 32 
<3-'32p = the thermal neutron cross section for the reaction P(n,^) P 

N32p = the number of atoms T present 
The ratio of both fluxes equals now 

'^f = "32scr32p N32p ^ Act32g<-'-32p N32p 

f D31p'~32s N32S Act32p0~32s ^325 
th 

32 
where Act32c, = the relative count rate of 'P produced by (n,p) reaction on sulphur 

Act32p = the relative count rate of P produced by (n,'/) reaction on phosphorus 

Experimental: Simultaneously known amounts of analytical grade (NH,)„SO, 

and (NH,)„PO, were irradiated at both irradiation sites used (pneumatic tube and 

reactor pool). Seven days after irradiation the salts were dissolved in diluted 

nitric acid and adequate amounts were counted by liquid scintillation using the 

V 32 

Cerenkov effect. P is a pure beta emitter with a maximum beta energy of 1.71 Mev 

35 
and a half life of 14.2 days. The interfering activity due to S produced by 

34 
(n,jj^) reaction on S is a pure beta emitter with a maximum beta energy of only 

•/ 

105 Kev and a half life of 85 days . This energy is too low to generate Cerenkov 

scintillation. To check the purity of the samples a decay curve analysis was 
performed for 3 months. For all the samples half- lives ranging from 14.0 to 14.4 
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days were found and no indication of the presence of shorter or longer lived isotopes 

was detected. The results shown in Table 9, are calculated using the following 

31 32 
generally accepted values for the cross sections (14,15). P(n,p) ^cr'=0.19 bam 

^■^S(n,Y) ^\^= 0.005 bam 



Irradiation 
Site 


Activity 

mg P 
c/min 


Activity 

mg S 
c/min 


Ratio 
Activity (£)_-.- 

Activities) 


Neutron 
Flux Ratio 
Thermal /Fission 


Pneum. Tube 
Reactor Pool 


206,480 
4,060,200 

, ~~~ _™_ 


8,823 
295,354 


23.40 
13.75 


7.6 
4.5 



Table 9. Determination of thermal to fission flux ratio. 

Calculation of the interferences: Based upon the experimentally found ratio 
of the fission to thermal neutron- fluxes and on the mean neutron fission cross 
sections as compiled by Neuerht and Zyp correction factors for possible interfer;- 
ences were calculated (Table 10- column 5). This correction factor equals the 
weight fraction of the interfering element which has to be subtracted from the 
element determined. Thus the corrected amount of element M present is given as 



^=^- ^E Vl. 



whe 



re M - the amount of element M experimentally found 

I = the amount of interfering element I experimentally found 

E 



Sl.I. 



= the correction factor for interference of element I on the deter- 



mination of element M 
The concentration ratios of interfering to interfered elements which give 
rise to an interference of 1% were also calculated. These ratios may be compared 
to concentration ratios found in actual air samples (Table 10, columns 6 and 7). 
It is obvious that in these samples a large correction 05%) can only be expected 
in the analysis of magnesium, which is interfered by a threshold reaction on 
aluminum. Therefore this low interference was also determined experimentally by 
irradiating a solution of pure A1(N03)3. It appeared that 100 micrograms of aluminum 



-29- 



Element 
Determined 


Interfering 

Element 

I 


Reaction 


Irradiation 
Site 


Correction 
Factor 

Sl.I. 


Concentration 
ratio 
resulting 
in 1% inter- 
ference I/M 


Actual 
Concentration 
ratios in air 
I/M 


Na 


Al 


Al (n , of ) Na 


Pneum Tube 


0.00015 


65 




2-8 


Na 


Mg 


24„ / n24„ 
Mg(n,p) Na 




0.00023 


43 


1-10 


Mg 


Al 


27.,, .27„ 
Al(n,p) Mg 


II 


0.15 


0.067 


0.5-5 


Mg 


Si 


30-., .27^ 
Sx(n,*<^) Mg 




0.00002 


490 


1-10 


Al 


Si 


2Ssi(n,p)2^Al 




0.002 


4.9 


0.2-4 


Al 


P 


31p(n,^)2«Al 


" 


0.00055 


18 




K 


Ca 


^'ca(n,p)^\ 


React Pool 


0.0013 


7.6 


0.5-5 


Sc 


Ti 


^^i(n,p)^^Sc 




0.000008 


1,250 


150-1500 


Cr 


Fe 


54^ , .51_ 
Fe(n,of) Cr 




0.000014 


720 


50-500 


Mn 


Fe 


56„ , n56,^ 
Fe(n,p) Mn 


Pneum Tube 


0.0000083 


1,200 


10-50 


Mn 


Co 


^^Co(n,c:<)^^Mh 




0.0000014 


7,200 


0.005-0.05 


V 


Cr 


52„ , .52„ 
Cr(n,p) V 




0.00002 


510 


0.5-5 



Table 10. Nuclear interferences due to threshold reactions. 



generates a false amount of magnesium equal to 20 + 8 microgram. This results 
in a correction factor K^ ...= 0.20 + 0.08 which is in fairly good agreement 
with the calculated value of K^ _... = 0.15. Magnesium concentrations are there- 
fore corrected with this factor (0.20) in actual analysis. 
7. Neutron Flux Mapping 

In the pool of the reactor up to 15 samples can be irradiated simultaneously 
in a polyethylene bottle with a diameter of 4 cm and a height of 8.5 cm. For the 
calculation of the concentrations it is assumed that all the samples receive the 
same neutron flux as the standard sample. However in the immediate neighborhood 
of a finite reactor core an important vertical and horizontal neutron flux gradient 
may exist. The samples are t hus liable to undergo a different flux. 

Experiment: An experiment was devised to determine the importance of this 
effect at the irradiation site used. An aluminum wire was stretched along the 
walls of the polyethylene bottle. After irradiation it was cut into pieces of 1 
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to 2 cm. (see Figure 9). These pieces were carefully weighed, counted and the 

specific activities were related to the neutron flux at the respective sites. 

72 71 72 

The Ga activity produced by the reaction Ga(n, 5') Ga on the gallium impurity 

24 
was counted as a measure of the thermal neutron flux and the Na activity produced 

27 24 
by the threshold reaction Al(n,t>s.) Na on the aluminum matrix was counted as 

a measure of the fission flux. These isotopes were counted by selecting their 

most prominent photopeaks respectively at 834 and 1368 Kev by means of a single 

channel analyzer coupled to a Nal detector. The result of one of these experiments 

is illustrated in Figure 9. As expected the fast neutron flux is more subject 

to gradients. Indeed a further thermalization sums up with the geometrical 

effect. The total horizontal gradient (4 cm) appears to be between 6 to 10% for 

the thermal and 25 to 35% for the fast fission flux. The total vertical gradient 

(8 cm) on the other hand is between 16 and 21% for the thermal and fission flux. 

All the elements, except one (Ni) are determined from mainly thermal neutron 

produced isotopes. However also the thermal gradients may cause non-neglectible 

errors. Therefore, only the lower 4 cm of the bottle are filled with one row of 

vertically stacked samples, with the standard in the middle. This largely reduces 

errors due to vertical gradients. Errors due to horizontal gradients are minimized 

by turning the bottle over approximately 180° at half of the irradiation time. 
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ll I. Computerized Data Reduction 
1. Introduction 

From the foregoing discussion it can be seen that complete analysis of Ge(Li) 
gamma spectra involves the steps indicated in Fig. 10. Data from the storage 
volvmie must be transmitted to the appropriate processing device. In the case of 
hand analysis this would involve a paper tape listing of the spectrum channel by 
channel or an oscilloscope display. Paper tape printout of 4096 channels takes 
minutes even with high speed printout devices, while writing the data on magnetic 
tape for later processing requires about 30 seconds or less. For computer process- 
ing, spectral data stored magnetic tape must be read into fast memory core storage. 

Next, all statistically significant peaks in the spectrimi must be located. 
Doublets should be recognized and Compton edges should be discriminated against. 
While this process of peak recognition may be easily accomplished by eye, it is 
a non-trivial problem for computer. Various computer techniques have been devised 

/-I (:\ 

such as the second difference method of Mariscotti, , cross-correlation analysis 

/■-I yN /-| Q\ 

of Black, the tangent method of Gunnink and others. In general, while 

all methods are successful in locating large isolated peaks, they may fail to work 

reliably in resolving doublets, discriminating between real peaks and Compton 

edges and in identifying small peaks which are marginally above the Compton continuum 

baseline. 

A problem of comparable difficulty is that of determining net peak areas. 
Analysis of peaks which are large in comparison with a relatively flat baseline 
can be accomplished by summation over a specified interval centered on the peak 
maximum with subtraction of baseline counts estimated from channels on each side 
of the peak as mentioned in Chapter 2. In the case of peaks located near Compton 
edges where the continuinn changes rapidly from channel to channel, this procedure 
can still be used successfully provided a judicious selection of channels is made 
for the estimate of the continuum under the peak. Difficulties arise in programming 
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(19) 
computers to make the appropriate choices. Ralston and Wilcox have developed 

a procedure for calculating the baseline using an iterative spectrum smoothing 
routine discussed below. This baseline spectrum is probably a better approximation 
generally to the Compton continuum than the so-called straight line approach above 
but is limited in its ability to handle cases of broad peaks located near Compton 
edges. Additional problems arise in estimating net peak areas from resolved 
doublets. If the peaks are of comparable height, areas of individual components 
may be proportioned on the basis of height with limited accuracy. A more satis- 
factory approach corresponds to the method of Covell developed originally for 
Nal spectra in which peaks are fitted by Gaussian functions to obtain resolved 
areas. If, however, a small peak resides on the shoulder of a larger one, the 
calculated net peak area of the smaller member is sensitive to the functional rep- 
resentation of the larger peak in the vicinity of the shoulder. At this laboratory 
a skewed Gaussian function combined with a linear background term has been used 
with some success in this type of situation but the increased running time in 
routine analysis makes it desirable to use alternative isolated peaks whenever 
possible. Fortunately, in Ge(Li) spectra, overlapping peaks are the exception 
rather than the rule so that alternative peaks can usually be found. 

The correspondence between peak position and gamma energy, typically expressed 
in terms of an N order polynomial, enables the computer to locate peak positions 
if gamma energies are supplied or to identify radioisotopes possibly contributing 
counts under a given peak by various detection processes. From a knowledge of 
the experimental conditions it is generally obvious which isotope is primarily 
responsible for the presence of a given peak, but a thorough search of a gamma 
energy-isotope library can be helpful in identifying possible interferences which 
might otherwise be overlooked. In routine analyses of spectra with the same set 
of peaks recurring, it is of course unnecessary to use a peak Identification option 
with each spectrum. 
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Net peak areas must be converted to isotope weights by means of appropriate 
factors which include a correction for the reduction in the ninnber of counts due 
to the dead time of the multichannel analyzer. Most modern multichannel analyzers 

have "live time" clocks which correct for Instrumental dead time by extending real 

(21) 
time counting periods. Gavron has pointed out that in addition to the fact 

that these clocks are not always accurate, they can not make the proper dead-time 

compensation for short lived samples where no simple relation exists between the 

(21) 
real time and the fractional dead time. An algorithm has been developed to 

handle this problem but it is desirable to keep counting rates sufficiently low 
that corrections having only approximate validity are unnecessary. In addition to 
conversion of net peak areas, the computer should make appropriate cooling and 
counting time decay corrections, fliox normalizations, subtractions for interfer- 
ences as indicated in Ch. 2 and calculate confidence levels in reported values. 

It is readily seen that hand extraction of the useful information content 
of gamma spectra is a tedious and time consuming endeavor subject to a variety 
of human errors. When spectral data are complex or when a large number of specjtra 
must be analyzed, as in the case of environmental sampling, it is especially de- 
sirable that data reduction be accomplished efficiently via computer processing 
with minimal human intervention. Because of the complexity of gamma spectra, 
computer programs have not yet been successful in taking over the entire job of 
data reduction, but the greatest share of the work can in fact be done automati- 
cally. 

Numerous computer-based spectrum analysis methods are reported in the liter- 
ature, especially for use with Nal spectra. The FORTRAN IV programs developed 
here for analysis of Ge(Li) spectra obtained from radionuclide mixtures present 
in suites of irradiated environmental samples, provide a broad data reduction 
capability. The operations outlined in Fig,, 10 are written in terms of subroutines 
which may be incorporated into main programs meeting the specialized requirements 
of the investigator. 
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In the following section, details connected with storage of spectral data on 
magnetic tape and its subsequent retrieval by the computer are considered. Suffi- 
cient procedural information is presented to make this section hopefully a useful 
guide on how to get data from the analyzer into the University of Michigan IBM 
360/67 computer. In section 3 the available subprograms and short main programs 
are described while in the final section the main program in current use with the 
analytical procedure outlined in Chapter 2 is described in detail. 
2. Tape analysis: 

A. Tape structure . Each spectrum is written on tape either as a set of 
N/256 records each of which is 258X6 BDC characters long or as one continuous 
record (N+2)X6 characters long where N=1024, 2048, 4096 channels. The writing 
mode is determined by a switch labeled (IRG/256) in the back of the ND 2200 series 
magnetic tape module. Existing programs expect the tape to be formatted in the 
latter way as one large record per spectrum. In this format the tagword is written 
twice at the beginning of the long record, while in the blocked format it is written 
twice at the beginning of each segment of 256 channels as shown in the ND 2200 
series magnetic tape manual. For each channel, counts are represented on tape 

by a string of six BCD characters. Since tagword values have at most four digits, 
two of the six characters for the tagword contain no useful information. 

B. Tape submission . Tapes should be submitted at the batch input window at 
the University of Michigan Computing Center, North University Building, Main 
Campus. Two labels are affixed physically to each reel submitted: (1) Computing 
Center ID(CCID) of the form Gxxx (e.g. G123) and (2) users tape ID (UTID) concocted 
for esoteric reasons by the user (e.g. SOURCE). In exchange for a tape the user 
receives a green card receipt having the appropriate Gxxx punched in the first 
columns. So as not to forget the association between the CCID and UTID it is 
desirable to write the UTID on the receipt card. Tapes so entered in the C.C. 
users library may be retrieved at the output window at a later time in exchange 
for the green card. Tapes may stay in the library so long as they are used at 
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least once a month. According to Computing Center personnel, tapes may not be 
retrievable if they are not used for several months or if the green card is lost! 
The UTID may be permanently associated with the tape while a new CCID will be 
assigned if the tape is withdrawn and subsequently resubmitted. 

The green card itself is never used as part of a series of program input 
cards. The information about v*iat tape to use is supplied on a $RUN *MOUNT 
command card described below. For the protection of the user, the CCID and UTID 
are intercompared before the tape is mounted by the operator. If the association 
is incorrect, the tape will not be mounted. 

C. Tape Inventory . Provided there are no parity or IRG "noise" errors on 
a tape, the string of BCD characters representing a spectrum can be translated 
into bytes and correctly grouped into Integers giving counts per channel by the 
usual formatted FORTRAN read statement. Often the original tape has errors and 
very short records sometimes arising from user intervention in the write operation 
at the wrong moment. The simple read statement will not handle such troubles 
for, upon detection of errors, control is passed back to MTS. 

A possible solution is to prepare an error-free tape using the MTS routine 
*TAPECOPY (Computing Center memo M 102) which can copy all records including 
those with parity errors. For faulty records the parity is re-calculated so as 
to be correct on the new tape. This program will not give information about the 
tagword and record length associated with each record. 

The FORTRAN tape inventory program INVEN (p. 94) does have this capability 
but will reject any record which Is not of correct parity. In fact the program 
as presently constituted will transcribe only records which are of correct parity, 
are "noise" free and of correct length (I.e. a multiple of 1024 channels). These 
more stringent requirements are adopted because it is not always clear how the 
record Is in error if the parity is wrong of if the record is of Improper size. 
An important advantage of having a tape Inventory, Is that it allows one to search 
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for a desired record on the basis of its position among many records rather than 
on the basis of its tagword. Even though it is desirable and customary to write 
record sequences with automatic tagword incrementing, it sometimes happens because 
of the human touch, that more than one record will have the same tagword, in which 
case record order is of clear importance. 

The input card sequence for INVEN is shown on page 38. 

Card 3: Beginning of MOUNT instructions. 

Card 4: Mount tape with CCID G123 on a seven track unit (there are only two 
at the Computing Center), assign pseudo device name *R* to it and convert (N) , 200 
char, per inch (2), even parity (E) characters to bytes using a tape buffer 
storage size of 25,000 characters. The UTID is SOURCE. By default (i.e. omitting 
the statement ring = out) the tape ring is taken out which means no writing can 
be done on this tape (which is the original after all!) The Computing Center 
operator will put in or remove rings as per instruction on the $RUN *MOUNT card. 
For additional information on *M0UNT consult Computing Center memo M67, M68 and 
CC News #145. ^^^-^ 

Card 5: Same as 3 but one writes on this tape. Hence ring = IN. UTID = CLTP6 

Card 6: End of MOUNT instructions. 

Card 7: Load and run the compiled main program +MVL in the file INVEN. 
FORTRAN uses integers instead of pseudo device names; hence the correspondence, 
2=*R*, 4=*W*, is made at load time. The program reads data from unit 2 and writes 
on tape mounted on unit 4. 

Card 8: Name of the tape, perhaps UTID, written on new tape inventory sheet. 
(2A4) format. 

Card 9 : Any additional information to be printed out on the new tape inven- 
tory sheet (20A4 format). 

Card 10: N, NSK in 2I4 format. 

The inventory begins after NSK records on the source tape have been skipped 
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and continues transcribing through record with tagword N before quitting. In the 
automatic incrementing mode this number can usually be equal to the last tagword 
written on the tape. If automatic incrementing is not used, a dummy spectrum 
should be written at the very end with N=9999. 

Program Notes: Details can be understood by inspection of the FORTRAN listing 
(p. 94). A unformatted read subroutine is used in which the modifier m is defined 
(hex: zAOOOOOOOO) so as to transmit error messages to the program rather than to 
MTS. A detected error will result in repeated read-back-space operations by the 
tape drive. This behavior is curtailed by invoking a one record SKIP operation 
(Computing Center memo M80) whenever more than four unacceptable files (too short 
usually) , four end-of-file marks or two peinnanent errors are detected in succession. 
This apparently arbitrary approach has greatly decreased the inventory time with- 
out impairing the efficiency of record retrieval. The most frequently encountered 
error is parity, code = 16. In the case of a tape drive failure, code = 24, the 
inventory will terminate. Experience has shown that if this error occurs, one 
should merely resubmit the job without changing any cards. Usually the second 
time around this error, which may arise because of the marginal compatibility of 
the Kennedy and the 360 tape drives, does not occur. Fortunately the code 24 error 
rarely occurs. 

Having produced a new tape, the progremi rewinds it and inspects it in 
essentially the same way for errors. Thus an error-free inventory of this tape 
indicates that all is well for subsequent data analysis. Cost of the inventory: 
100/4096 channels. 

3. Output: Part of a source inventory is shown on page 40. EOF=End of file, 
FLN-File number, LEN=number of characters in the record, CHNLS= corresponding 
number of channels, CMRD=cumulative number of read operations, CRN=index of 
acceptable records written on new tape. Of the 42 records on this tape only one 
(Tg=0039) was found to be unacceptable. Thus the inventory of the new tape also 
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Input cards for Running the 
compiled Fortran Program INVEN. 



1. $SIGNON SN87 T=IM 'NAME' 

2. PASSWORD 

3. $RUN *M01INT 

4. G123 ON 7TP *R* M0DE=2EN SIZE=25000 'SOURCE' 

5. G308 ON 7TP *W* M0DE-2EN RING-IN SIZE=2500 'CLTP6' 

6 . $ENDFILE 

7. $RtJN INVEN 2=*R* 4=*W* 

8. CLTP6 

9. WRITTEN ON 9/26/69 

10. 00420000 

11. $S1G 
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listed on p. 40 displays the tagword and record length for 41 records. All 
extra end-of-file marks between records have been deleted. 

Shown on p. 41 is part of a source inventory for an unusually bad case. The 
program detected two parity errors in a row after Tg=0041, decided thus to skip 
one record, immediately found a record far too short followed by another record 
with incorrect parity. After skipping a region with too many end of file returns, 
it detected a series of extremely short records which in the case of "Tg"=0213 
it attempted to read 3 times before invoking skip. Notice that despite patches 
of trouble, no useful information was lost. The records were written in the auto- 
matic incrementing mode and each tagword is in fact present on the transcribed 
tape (CRN=14). Occasionally a desired spectrum is lost, as in the previous example 
(Tg 0039). The ratio of lost to saved spectra is about 1/50. It is desirable 
to write very important spectra twice to insure their retrievability. The 
Kennedy tape recorder accepts 1200 foot reels which can store about 100 4096 
channel spectra. Thus one reel will take about 50 duplicate spectra of this size. 
3. Data Analysis Programs. 

Large multipurpose programs are built out of collection of subroutines having 
specific tasks, thus permitting rapid assembly of programs to fit specialized 
needs. In this section the subroutines are described as well as several short 
main programs which make use of them. 

A. Tape read and list (RDLST) . A set of routines read the tape (either the 
new tape produced by the inventory program or an error free original tape) and 
print out the contents of each channel. The main program (p. 9$ (RDLST) makes 
use of five subroutines RDTAPE, MVL, WRTMT,, BASE, and WRTBS. As seen from the 
FORTRAN listing of RDTAPE (p. 95) the tape read format is 3A4, 32(12816). Reading 
in the first 12 characters as 4 byte words avoids complications connected with 
the special encoding of the tagwords. Namely: 
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Inventory Program Output 



COkT£MS OF TAPE BEIhG RFAO 



CONTEKTS Of NEULV HRlTFPh TaPE CLTP6 
hRITTEK OK 'i/lb/bt . 







FLN 




TGW 


COCl 


LEN 


12303 


CHNLS 


2040 


CMRO 


3 


CRN 


1 


FILE 


hO. 


1 


TGbC 


0001 


INCTH« 


1230C 


CHR 


CHNLS 


204B 


CUM 


ROS 


ECF 


3 


























FILE 


NO. 


2 


TGWC 


CC02 


LNGTM" 


12300 


CHR 


CHNLS 


2048 


CUM 


ROS 


ECf 


4 


























FILE 


hC. 


3 


TGkC 


0003 


LNCTH- 


12300 


CHR 


CHNLS 


2048 


CUH 


ROS 






FIN 




IGw 


CCC2 


lEN 


12 300 


CHNL S 


2048 


CHRfl 


6 


CRN 


2 


FILE 


hO. 


4 


TGhD 


0004 


LNGTH" 


12300 


XHR 


CHNLS 


2048 


CUM 


ROS 


ECF 


5 


























FILE 


NO. 


5 


TGWC 


O005 


LNGTM- 


12300 


CHR 


CHNLS 


2048 


CUM 


ROS 


ECF 


6 


























FILE 


^0. 


6 


TGkC 


0006 


LNGTH» 


12300 


CHR 


CHNLS 


2D4R 


CUH 


RDS 






FLN 




TGW 


OCO) 


LEN 


12300 


CHNtS 


2048 


CMRO 


9 


CRN 


3 


FILE 


NO. 


7 


TGWC 


0007 


LNGTH- 


1230C 


CHR 


CHNLS 


2048 


CUM 


RDS 


ECF 


7 


























FILE 


NC. 


e 


TGWC 


0008 


LNGTH- 


12300 


CMR 


CHNLS 


2048 


CUM 


RDS 


ECF 


a 


























FILE 


hC. 


<; 


TGkC 


0009 


LNGTM« 


12 300 


CHR 


CHNLS 


204 B 


CUM 


RDS 


ECF 




FIN 




TGW 


COOA 


LEN 


12300 


CHNLS 


204ft 


CM«D 


12 


CRN 


4 


FILE 


NO. 


u 


TGhC 


COlO 


LNGTM. 


12300 


CHR 


CHNLS 


2048 


CUH 


ROS 




























FUF 


NO. 


11 


TCkC 


0011 


LNGTM. 


12300 


CHR 


CHNLS 


204 B 


CUM 


ROS 


EC 


10 


























FILE 


NC. 


12 


TGhC 


0012 


LNGTH* 


12300 


CHR 


CHNLS 


2048 


CUM 


RDS 






FLN 




TGW 


CCC5 


IFN 


12 300 


CHNLS 


2048 


CHhO 


1& 


CHN 


5 


FILE 


NC. 


13 


TGWC 


0C13 


LNGTM- 


12300 


CHR 


CHNLS 


2048 


CUM 


ROS 


ECF 
tOF 


11 


























FILE 


KC. 


t« 


TGkC 


C014 


LNGTH- 


12300 


CHR 


CHNLS 


2048 


CUM 


ROS 


12 


























FILE 


NO. 


15 


TGWC 


0015 


LNGTH» 


12 30C 


CHR 


CHNLS 


2048 


CUH 


ROS 


ECF 

ECF 




FLN 




1GM 


C006 


LEN 


12300 


CHNLS 


2048 


CMRD 


IB 


CRN 


6 


FILE 


KC. 


16 


TGWD 


0016 


LNGTM- 


12300 


CMR 


CHNLS 


2048 


CUH 


ROS 


li 


























FILE 


KC. 


17 


TGhD 


OOIT 


LNGTH- 


12300 


CHR 


CHNLS 


2048 


CUM 


RDS 


U 


























FILE 


NO. 


18 


TGWO 


OOIB 


LNGTH- 


12300 


CHR 


CHNLS 


2048 


CUM 


RDS 


ECF 




FLK 




TGW 


0007 


LEN 


12300 


CHNLS 


2048 


CHRO 


21 


CRN 


7 


FtLC 


NC. 


IS 


TGkO 


0019 


LNGTH- 


12300 


CHR 


CHNLS 


2048 


CUH 


ROS 


15 


























FUF 


NO. 


20 


TGhD 


0020 


LNGTH- 


12300 


CHR 


CHNLS 


2048 


CUH 


RDS 


ECF 


16 


























FILE 


NO. 


21 


TGKC 


0021 


LNGTM" 


12300 


CHR 


CHNLS 


2048 


CUM 


RDS 


ECF 




FLN 




TGU 


COCH 


LEN 


12300 


CHNLS 


2048 


CHRO 


24 


CRN 


S 


FILE 


KC, 


22 


TGhO 


0022 


tNGTM« 


12300 


CHR 


CHNLS 


2048 


CUH 


RDS 


17 


























FILE 


SO. 


23 


TGWC 


0023 


LNGTM. 


1230C 


CHR 


CHNLS 


2048 


CUM 


ROS 


Ifi 


























FILE 


NO. 


24 


TGwn 


0024 


LNGTM- 


12300 


CHR 


CHNLS 


2048 


CUM 


ROS 






TLN 




rcw 


CCC9 


lEN 


12300 


CHNLS 


?04d 


CMRD 


27 


CRN 


9 


FILE 


NC.' 


25 


TCUC 


0025 


LNGTH- 


12300 


CMR 


CHNLS 


2048 


CUM 


ROS 


ECF 
ECF 


19 


























FILE 


NQ. 


26 


TGhC 


0026 


LNCTM- 


12300 


CHR 


CHNLS 


2048 


CUM 


ROS 


20 


























FILE 


NO. 


27 


TGWC 


0027 


LNGIM- 


12300 


CHR 


CHNLS 


2048 


CUH 


ROS 






FL^ 




TGH 


CCIC 


LEN 


12300 


CHNLS 


20Aa 


CMRD 


30 


CRN 


10 


FILE 


NO. 


28 


TGhr 


0028 


LNGTH- 


12300 


CHR 


CHNLS 


204B 


CUM 


RDS 


ECF 


21 


























FILE 


hO. 


It 


rr.wc 


0029 


LNGTH- 


12300 


CHR 


CHNLS 


2048 


CUM 


RDS 


22 


























FILE 


no. 


30 


IGhC 


00 30 


LNGTM- 


12300 


CMR 


CHNLS 


2048 


CUM 


ROS 






FLN 




TGK 


coil 


LEN 


12 300 


CHNLS 


20 4 S 


CMRO 


33 


CRN 




FILE 


SO. 


31 


TCWC 


0031 


LNGTH= 


12300 


CHR" 


CHNLS" 


2048 


CUM 


RDS 


tCF 


23 


























FILE 


NO. 


32 


TGkfl 


0O32 


LNGTH> 


12300 


CHR 


CHNLS 


2048 


CUM 


RDS 


Etf 


2* 


























FILE 


NO. 


3' 


TGhD 


0033 


LN'GTh- 


12300 


CHR 


CHNLS 


2048 


CUM 


ROS 


ECF 




FIN 




TGW 


C012 


LEN 


12300 


CHNLS 


2048 


CMRD 


36 


CRN 




FILE 


NC. 


34 


TGWC 


0034 


LNGTM» 


12300 


CHR 


CHNLS 


2048 


CUM 


RDS 


29 


























FILE 


hC. 


35 


TGWC 


0035 


LNGTH- 


12300 


CMR 


CHNLS 


2048 


CUH 


ROS 


CCf 


26 


























FILE 


NO. 


36 


TGhO 


0036 


LNGTH* 


12300 


CHR 


CHNLS 


2048 


CUM 


RDS 


ECF 
ECF 


27 
28 


FLV 




TGW 


ecu ■ 


LFN 


123 JO 


CHNLS 


204a 


CMRD 


39 


CRN 




FILE 


NC. 


37 


TGHC 


0037 


LNGTH. 


12300 


CHR 


CHNLS 


2048 


CUM 


ROS 


























„^FILE 

^•fue 


NC. 


38 


TGkC 


0038 


tNGTM» 


12300 


CHR 


CHNLS 


2048 


CUH 


ROS 


























NO. 


3<: 


TGWD 


0040 


LNGTHa 


12300 


CHfi 


CHNLS 


2048 


CUH 


RDS 


ECF 
ECF 




FLN 




TGW 


con 


LEN 


12300 


CHNLS 


2048 


CHRO 


42 


CRN 




FILE 


NO. 


AC 


TGhC 


0041 


LNGTM- 


12300 


CHR 


CHNLS 


2048 


CUH 


ROS 


29 
30 


























FILE 


NO. 


41 


TGkC 


0042 


LNGTH. 


12300 


CHR 


CHNLS 


2048 


CUH 


RDS 


ECF 


31 


FLN 




TGU 


C015 


LEN 


12300 


CHNLS 


2048 


CMRO 


45 


CRN 






















- 






ECF 

ECF 


12 

13 

14 


FIN 




TGW 


C016 


LEN 


12300 


CHNLS 


2048 


CMRD 


48 


CRN 




























CCF 


















































ECF 


l^ 


FLN 




TGW 


CCI7 


LEN 


12300 


CHNLS 


2048 


CMPD 


51 


CRN 




























ECF 


M 


FIN 




TGH 


C018 


LEN 


12300 


CHNLS 


204« 


CMRO 


S4 


CBN 


16 



























« \V,c 







FIN 


3? 


ECF 


79 






ECF 


to 










FIN 


40 


ECF 


81 






CCF 


82 










FLN 


41 


ECF 


83 






ECF 


84 










FLN 


42 



TGW CC39 LEN 1444 CHNLS ?39 CHAD 117 

TGW C040 LEN 12300 CHNLS 2043 CMRO 120 CRN 39 

TGW C041 LEN 12100 CHNLS 2049 CMRO 123 CRN 40 

TGW C04? LEN 12300 CHNLS 2048 CMRD 126 CBN 41 



v- siM»«r ^ecpM>! 
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A Source Tape Listing of an Unusually Bad Tape 







COMTFNTS 


OF TAPe BEING 


BEAD 


































































■ 


■OF 


1 


FLN 


1 


T6W 0035 


LEN 


24588 


CHNLS 


409i!) 


CMRD 


1 


CRN 


1 


EOF 


? 


FtN 


Z 


row 0036 


LEN 


24588 


CHNLS 


4096 


CMRD 


4 


CRN 


2 


;oF 

=0F 


3 
4 
























■OF 


5 


FLN 


3 


TGW 0037 


LEN 


24588 


CHNLS 


4096 


CMRD 


7 


CRN 


3 


ioF 


6 


FLN 


<t 


TGrf 0038 


LEN 


24588 


CHNLS 


4096 


CMRD 


10 


CRN 


4 


roF 

:0F 


7 
8 
























EOF 


9 


FLN 


5 


JGA 0039 


LEN 


24588 


CHNLS 


4096 


CMPO 


13 


CRN 


5 


;0F 
iOF 


10 
11 
























iOF 


12 


FLN 


6 


TGW 0040 


LEN 


24588 


CHNLS 


4096 


C«*RO 


17 


CRN 


6 



IOF 
IOF 


13 
14 






























IRR 


1 


CD 


16 


FLN 


1 


TGW 


0041 


LEN 


24588 


CHNLS 


4096 


CHRO 


21 


CRN 


7 


IRR 
fOF 


2 
15 


CD 


16 


SKIP 
























IOF 


16 






FLN 


■ 8 


TGW 


0042 


LEN 


317 


CHNLS 


50 


CMRD 


25 






EOF 


17 






FLN 


9 


TGW 


0042 


LEN 


24588 


CHNLS 


4095 


CMRD 


27 


CRN 


8 


tOF 
tOF 


18 
19 






























IRR 


^ 


CO 


16 


FLN 


10 


TGW 


0043 


LEN 


24588 


CHNLS 


4096 


CMRD 


31 


CRN 


9 


IRR 
OF 


4 

20 


■CO 


16 


SKIP 
























;0F 
EOF 


21 
22 






























[of 


25 


SKIP 




FLN 


11 


TGW 


0044 


LEN 


24588 


CHNLS 


4096 


CMRD 


38 


CRN 


10 


tOF 
EOF 


24 
25 






























fOF 
EOF 


26 
27 






FLN 


12 


TGW 


0045 


LEN 


24588 


CHNLS 


4096 


CMRD 


41 


CRN 


U 






FLN 


13 


TGW 


0046 


LEN 


24588 


CHNLS 


4096 


CMRD 


43 


CRN 


12 


50F 
;0F 


28 
29 






































FLN 
FLN 


14 
15 


TGW 
TGW 


004 7 
0002 


LEN 
LEN 


24588 
36 


CHNLS 
CHNLS 


4095 
4 


CMRO 
CMRD 


47 
48 


CRN 


13 










FLN 
FLN 


16 
17 


TGW 
TGW 


0213 
0213 


LEN 
LEN 


27 
27 


CHNLS 
CHNLS 


2 
2 


CMRO 
CMRD 


49 

50 






PUEER 


FILE SKIP 


FLN 


18 


TGW 


0213 


LEN 


27 


CHNLS 


2 


CMRO 


51 






EOF 


30 






FLN 


19 


TGW 


0048 


LEN 


81 


CHNLS 


11 


CMRO 


52 







FLN 20 TGW 0048 LEN 24588 CHNLS 4096 CMRO 54 CRN 14 
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Tagword #1 Tagword #2 




II 12 13 IT 

3A4 16 etc. 

(13 = 1234) thus is the correct tagword In A4 format. The short subroutine MVL 
(p. 95) converts 13 to I format so that tagword input information from cards is 
conveniently compared. Subroutine RDTAPE then will read one record of length 
NR+1 channels. NR '(e.g.: 1023_, 2042 o^ 2095^) must be supplied by the program 
while the tagword, 13 (in alphameric mode), the [ ^] time of data accumulation, 
IT, and the data array IR(J) (J=1,NR) are returned. IT is stored in the "zeroth" 
channel so that only NR channels contain pulse height data. Subroutine WRTMT 
(p. 95) will print out the values 13, IT and (IR(J) , J=l, NR) obtained previously. 
BASE and WRTBS are described below. 

The input card sequence required to print out portions of five of the forty- 
one spectra listed in the inventory of CLTI6 p. 44 are shown below. Six values 
entered on one card are needed for each spectrum listed; NTG, NREC, NSK, IS, IE 
and NBS in 6l4 format. NTG = the tagword sought, NREC = the total number of 
channels in the spectrum (i.e. 1024, 2048, or 4096), NSK = the number of records 
to be skipped starting with the present position of the tape, before searching 
for the target record (i.e. that with the desired tagword). The search operation 
consists in successive formatted reads and comparison of the NTG with the tagword 
read, 13. If the tagword sought is less than the tagword read, a mismatch is 
detected and the program will stop. This feature can easily be changed so that 
the tape will backspace an appropriate number of records in its searching but 
this would be desirable only if automatic incrementing has been completely success- 
ful. In this connection it is important to note that the operation of skipping 
records is much faster than the formatted read operation. Hence it is desirable 
to choose N3K so that the very next record is the target record. This can be 
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done on the basis of the tape inventory which gives the correspondence between 
the record index and tagword on the newly written tape. IS and IE are respectively 
the first and last channels of the spectrum to be printed out. The listing gives 
500 channels per page of output. NBS <_ for a simple listing. The case for NBS 
> is discussed below. 

The first card in the data set. Card 7,, causes the last 1048 channels of 
spectrum with tagword 1 to be printed out. Card 8 causes the first 2000 channels 
of the spectrum with tagword 7 to be printed. In this case 5 spectra are skipped 
prior to searching for 7. As mentioned previously if NSK had been taken = the 
same output would have been obtained but the running time would have been in- 
creased because spectra 2 through 6 would have been read in in order to find 7. 
A portion of the first page of output for a spectrum is shown below. The 
time is labeled "live-time" by the program in as much as this is the common oper- 
ational mode of the analyzer but the actual meaning of this number is of course 
determined by the position of the live- time switch on the analyzer itself. The 
effects of cards 9 and 10 can be seen by inspection of the inventory on p. 40. 
A blank card will cause RDLST to stop. 

B. Base spectrum generation and peak integration . The subroutine BASE 
(p. 95) provides a measure of the pulse height spectrum baseline on which peaks 
rest, thus enabling the computer to make a determination of net peak areas. The 
technique is described in detail elsewhere (Ralston Se Wilcox 19) but, briefly, 
the program generates from the original spectrum, S°, with n. counts in the ith 
channel, a smoothed spectrxim S*-" in which the number of counts in the ith channel, 
n., is taken as the average over an interval of width 2£+l centered about i. 
Thus i -iji-tA 

_ \ 




S'-i-A 
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Card inputs to RDLST 



1. $SIGNON SN87T=IM 

2. PASSWORD 

3. $RUN *MOUNT 

4. G308 ON 7TP *R* M0DE=2EN RING=OUT SIZE=25,000 'CLTP6' 

5 . $ENDFILE 

6. $RUN RDLST 4 = *R* 

7. 000120480000100020480000 

8. 000720480005000120000000 

9. 003820480032050010000000 

10. 004020480000000120000000 

11. 004220480001070009000030 

12. BLANK CARD 

13. $SIG 
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S° and S^ are compared. Whenever n, exceeds n. by a standard deviations (i.e. 
by a*/ n. ), n. is replaced by n, thus generating a new spectrum S-'-. 

The above steps of smoothing and replacement are repeated in turn using S 
in place of S° and so on for n cycles with ever-diminishing values of a. (<^-ii 
<a.). The resultant spectrum S^ is taken as the baseline. In subroutine BASE 
the initial value of a is taken as 9.2, decreasing in ten (i.e. n=10) equally 
spaced steps to 0.2. £ is taken to be 30 channels. The choice for £ results 
from a compromise between the use of a small £ (~10 channels) for which the 
Compton edges are well-followed and a large value which gives a better represen- 
tation of the baseline under unusually broad peaks. An example of the effects 
dependent on the choice of £ is shown in Fig. 11. The spectrum was taken at a 
very high counting rate in order to exaggerate peak widths. With £=30 channels, 
the calculated baseline is largely insensitive to peak width variations encountered 
in actual measurements. 

In the plot of a typical Ge(Li) gamma Spectrum (Fig. 12). The smooth curve 
follows the upper limit of the calculated baseline. Portions of the spectrum 
which do not exceed 0.2 n. are not replaced in forming the baseline. Thus only 
under the peaks is the baseline smooth. Elsewhere it often has the value of the 
original data. The first and last (2 £+1) channels (61 in this case) are lost 
in the analysis. It will be seen that the baseline follows well the curve that 
would be drawn by eye except in the vicinity of Compton edges. While this program 
provides a means for doing routine analysis of peaks known to be away from Compton 
edges, it is the responsibility of the user to insure that this is in fact the 
case for his particular experimental conditions. Also note that BASE tends to 
underestimate the "baseline" associated with small peaks on the shoulder of a 
large one as in the case of peak labeled 10 situated on the shoulder of peak 11. 
Arguments of subroutine BASE (p. 95): 
L = sum half width = 30 
1ST = starting channel (usually = 1) 
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IND = final channel (usually - NR) 
IR = data array 

BS = returned base spectrum array 
The subroutine PKBS listed on p. 95, "integrates" the peak simply by summing 
counts over 2*NFWfl channels centered about channel ICH. Base spectrvmi counts 
obtained previously are sximmed over the s^mie interval. The subroutine returns 
SBS= the net counts = 




ll^u^-?sa)^ = N 



B 



the standard deviation, SD = / N + B, and the ratio = SD/SBS. In the expression 
for SD it is assumed that B has the Poisson- statistical uncertainty one would 
associate with real counts. For real peaks of about 3KeV FWHM, when the energy 
calibration is IKeV/ channel, the ratio SD/SBS as a function of NFW has a minimum 
around NFW=2. In order to compensate for possible line shape variations, NFW 
is typically taken = 3 which slightly increases SD in general (a few percent in- 
crease). PKBS requires the "peak" channel, ICH, as input. Thus any drift com- 
pensation must be done before calling PKBS. 

A subroutine WRTBS (p. 95) operates like WRTMT but lists the values of the 
base spectrum for each channel as well as the original data. In the read-list 
program RDLST, if NBS >0, a base spectrum will be calculated from channel IS, to 
IE with L=NBS (usually = 30) and WRTBS rather than WRTMT will be called. A sample 
portion of WRTBS output is shown on (p. 49) for the last spectrum on CLTP6. This 
spectrum is produced from card 11 in the input (p. 44) to RDLST. This print out 
was originally obtained in order to examine the behavior of the computed baseline 
in the vicinity of the 846.9 KeV Mn peak. It will be seen that the actual counts 
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In a given channel (e.g. the "peak" channel, 842) is the upper number of the pair 
(2144) while the lower one is the baseline count (261). Arguments of subroutine 
WRTBS: 

IS = starting channel (usually 1) 

IE = final channel (usually NR) 

BS = previously calculated base spectrum 

C. Plotting . A makeshift program has been written to product plots like that 
shown in Fig. 4. Limited to six consecutive plots, the program requires 5 values 
for each spectrum: NREC, ITG, NCUT, NSK, NCR. 

NREC = no. of channels 
ITG = tagword desired 

NSK = skip NSK records before looking for ITG 
NCR = do not calculate base 

>0 calculate base spectrum and plot 
NCUT = no. of channels to be plotted 
NREC = stop 
The channel axis must be <10 inches. For 4096 channels this gives FT = .0025. 

While the log counts axis can be up to 24 inches. But for a fit on 8 1/2" x 11" 

4 
paper, assuming 10 cts maximum counts in any one channel, FT = 2 gives a maximum 

height or 8" = 2 log lO'^. 

The plot shown in Fig. 4 was produced from the series of input cards shown on 

(23) 
p. 50. For a description of the MTS control cards, obtain the MTS plotting manual 

D. Peak search routine . 

A peak finding routine, PEAK, (p. 97), is available which has been tested and 
found to work well except in the case of doublets. This subroutine which in turn 
calls a series of subroutines FKZ, CIJWZ, SFWZ, locates peaks using a second dif- 
ference operator function described in detail by Mariscotti.^ "^ A FORTRAN listing 
of a test main program, PKTST, and the above subroutines is on p. 97. The twenty-five 



-49- 



o 

ro 
II 
w 

PQ 

Z 



o 

•rl 
4-1 

& 

o 



u 

Ci. 

w 

CO 

+-1 



H 
00 



e 

o 
u 

■M 
C 



4J 
M 
cd 






0< CO 
GO ^ 



oo 
u> in 

(M CM 



CO CO 



CO CO 
(M <M 



O — 

>-l 00 



•^ CO 



CO H 
O CO 



CO -t 

a 00 



0> » 



Psl (Nl 



(M CO 

ot- 



M (M 

ir\ in 



<0 O 

IM CVi 



CO 0> 

CM pg 



« <0 

CM CM 



r- f* 

N CM 



N CM 



f- 50 
O 1^ 
ft CM 



in CO 
o< N 

N CM 



h- i-i 
OCO 



ON <J> 
CM CM 



tM <-> 
f4 00 

m (M 



CO 00 
M CM 



-< O 

O CO 

rn (M 



n in 

m in 

IM CM 



!^ o 

oo 00 
IM (M 



5> 0> 
N CM 



1J> 0> 
IM CM 



OO 00 

tM CM 



o 
o 

00 




o 

00 




o 

<M 
CO 


o o 

00 00 
CM CM 




0> 00 
IM CM 




OO 
CO CO 
CM IM 


(M CM 




in in 
•o -c 

CM CM 




CM -I 
O CO 

CO «M 


O 00 
CO CM 




in in 

CM (M 




-* IIM 

m IN 


<»• CM 

•H OO 
tr CM 




o o 

(M CM 




IM 0« 
CM 00 
N»- IM 


<M CM 




CM CM 
CMCM 




O ID 

Nt IM 


O O 
00 00 
CM CM 




ON 00 

<M CM 




o m 
in IB 
CO IM 


in in 

CO 00 
CM IM 




0> 00 
IM (M 




O c*l 
O IB 
CO IM 


<M CM 




CO CO 
O 00 
m CM 




mm 
00 <s 
CM (M 


■O CM 
O 00 

m CM 




CM CM 




N IM 


00 CM 
<M CO 
<*l CM 




cM rsi 




o> in 

CM IM 


•-4 




«"4 

o 

00 




t-l 
CO 



IM 
CO CM 



m CM 
00 ^- 

fM IM 



<t CO 
CM CM 



CO ^ 
o r- 

(O CM 



» m 
CO h- 

CM CM 



nO <o 
o f- 

CO IM 



Of» 

or~ 

CO CM 



lO CM 



o o> 

lO IM 



m ON 
o r^ 

CO CM 



O CO 

r» <o 

<B CM 



CM •* 
ON ^ 
* CM 



m m 

CO nO 
* CM 



CM m 

nO « 

CO CM 



o <0 

CO CM 



0> .0 

IM CM 



m 00 
o nO 
lO CM 



CO ON 

" nO 

m CM 



N 0> 
CO o 
CO IM 



r^ o 

CO f» 
CO IM 



n 

c 



IM CM 
rg CM 
IM CM 



IM CM 
CM IM 



fH CO 

m 

IM CM 



o m 

V m 

CO IM 



nO 

h- m 

CM 



no m 

CM CM 



<t m 

m (M 



CO O 
<* nO 

m CM 



■4- nO 
>4CM 
CM 



•4' CM 

00 IM 



ON ON 

im4 w^ 
CM CM 



m m 

CM CM 
CM IM 



o o 

IM IM 
IM CM 



h- IM 
IM IM 



O O 
IM IM 



IM IM 
>*■ >* 
IM IM 



(O CO 
IM IM 



IM CM 
CM CM 

CM CM 



•O nO 
lO CO 
IM CM 



CM CM 
IM CM 



ON e> 

■-4 •-4 
CM CM 



ft f-> 
IM IM 



nO nC 
CM IM 
IM IM 



O O 
CO CO 
CNJ IM 



0> IM 
•*■ CO 
IM CM 



lO CO 
IM CM 



r-Nj- 
•* CO 
CM CM 



0> ON 
IM IM 
CM IM 



IM CM 
IM CM 



O O 
IM IM 



o 

00 
CO 


nO <0 

o o 

IM CM 


CO O 
<*• IM 

CM CM 



m m 
o o 

IM IM 



O -4 

m (M 

CM (M 



oo 

O O 
IM (M 



Nt Nt 

IM CM 
CM IM 



IM CM 
CM CM 



O O 
CM CM 



lO OO 
O CM 
IM CM 



0> f- 

m CM 

CM CM 



-50- 



Input Card Listing for Plot Program. 

$SIGNON SN87 T=IM 

PASSWORD 

$CEEATE PLOT SIZE=1000 

$RUN *MOUNT 

G308 ON 7TP *W* M0DE=2EN RING=OUT SIZE=25G00 'CLTP6" 

$ENDFILE 

$RUN *FORTRAK 



FORTRAN LISTING P. 96 



$ENDFILE 

$RUN-LOAD#+*SOURCE*+*PLOTSYS 4=*W* 9=PL0T 



COMPILED SUBROUTINES (RDTAPE, MVL BASE) 



$ENDFILE 

4096000818000006000 

BLANK CARD 

$RUN *CCQUEUE PAR=PLOT 

$RUN ^PERMIT PAR=PLOT 

$SIG 
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labeled lines shown in the plot (Fig. 12) were identified as peaks by the program. 
It will be seen that even a small peak on the shoulder of a large one is currently 
detected as in the case of pair 3-A or 10 and 11. Also Compton edges are recog- 
nized and excluded on the basis of behavior of the second difference function. 
The sensitivity of peak detection is adjustable. In the present case, peaks 
labeled (A) and (B) were not identified but would have been if the sensitivity 
had been increased. 

The peak finding routines have not as yet been incorporated into composite 
programs for routine use because in quantitative analysis of environmental samples 
the same set of peaks generally recur in each spectrum. Further information on 
this program is available on request. 

E. Component Identification . The subroutine IDCOM (p. 98) makes use of a 

(12) 
list of 243 gamma-emitting nuclides and their energies to determine which nuclides 

may contribute by various detection processes to a gamma spectrum peak observed 

in a channel corresponding to the energy E,. If the fwhm of the peak is 2*6E, a 

gamma ray of energy , Ey , will contribute provided 

E^ = E + 6E 
for complete energy absorption in the detector ("main" peaks) or 

Ej, - SllKeV = E + 6E when E >1022 KeV 
for single-escape peaks or 

Ey - 1022 KeV = E + 5E E >1022 KeV 
for double-escape peaks, or provided two gcimmas from the same nuclide sum to the 
correct energy 

^Vl "•■ ^>'2 " ^ - ^'^' 
The nuclide data is contained in a private file GAMLIB but is available on 

cards. It includes a complete list of slow-neutron produced (n,}*) nuclides as 

137 
well as some fast-neutron produced species and background nuclides such as Cs 

and K . Actually the file contains two lists : (1) A main list with nuclides 
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index, name, half-life (min.) and up to five energies of the most intense gammas 
emitted (primary gammas) and (2) an auxiliary list containing up to 100 secondary 

gamma energies. The division into primary and secondary groups is somewhat quali- 

(12) 
tative following the classification of Adams ani Dams. 

IDCOM requires a set of measured gamma energies, (ES, in KeV) the total 

number in the set, IPK, and the time (in minutes, TIME) elapsed between the end 

of irradiation and the start of counting. The subroutine expects to treat an 

entire spectrum. Thus if one finds no "main" peaks associated with the primary 

gammas of a certain nuclide, one certainly should not expect to find any S.E. , 

D.E., sum or secondary peaks from that nuclide in the spectrum. This might not 

be true if certain peak energies were missing from the set. IDCOM first makes 

a list of those nuclides for which there is in each case at least one main peak 

found and for which the elapsed time is less than ten half -lives. Other nuclides 

are not considered further. For each nuclide on this list the subroutine calculates 

S.E. , D.E. and sum energies from the primary gamma energies, obtains secondary 

peak energies from the auxiliary GAMLIB list and looks for equalities between 

these values and the energies ES + 6E. Each nuclide name and coded detection 

process attributed to a given peak is indexed and stored for eventual transmission 

out of the subroutine. Finally for all j with ES(j) <100 KeV, the subroutine looks 

1 ES ( •■''■^ 

for an ES(J ) satisfying the relation, 3 = 1.132, which holds for the Ka 

ES(j) 

and K„ X-ray pair. If such a pair is found, the Z of the associated atomic nucleus 

p 

is given by 

Z = f(E) = 28.553 + 0.9485 E - 0.003085E^ 

FWHM 
The value of the energy spread, 6Es — r— , is obtained from a function subroutine 

EF(E). IDCOM Arguments: (values returned by IDCOM; not intermediate values.): 

N(j) = number of nuclide-process attributes associated with the jth peak of energy 

ES(j), NI(K,J) = process code associated with the Kth attribute and Jth peak. Thus: 
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NI = 1 main peak 

2 double escape peak 

3 single escape peak 

4 sum peak 

5 secondary peak 

6 X-ray pair member 

IDX(K,J) = Nuclide index 1 £ IDX<243) associated with the K*^ attribute and J*^ 

peak. If NI = 6, IDX is not the isotope index, but rather the Z of the X-ray 

transition nucleus. IDX = 0.5 + f(E). 

NS(L, Loc) = Nuclide symbol for the L nuclide on the nuclide list. Two locations 

required for the symbol. L0C=1,2. 

NPP = total number of nuclides in the list. 

IP(L) = Nuclide index associated with the L nuclide. This value is needed for 

finding the correct secondary gamma energy list. 

INDX(L) = Actual number of primary gammas (excluding the 511 KeV gairana) associated 

■f* Vi 

with the L nuclide. 

NPKL (L) = Number of primary gammas (excluding 511 KeV) associated with the L 

nuclide for which a main peak was present. NPKL(L) <_ INDX(L). 

ER(L,k) = The primary gamma energies associated with the L nuclide. K=l, 5. 

If N(j) = the j peak is said to be "not attributed". 

Subroutine WRTID (FORTRAN p. 98). This subroutine will output information from 

IDCOM. It must be called for each peak, of energy ES(j). From a data array IZ, 

WRTID picks the correct atomic element symbol to go with the calculated Z when 

NI=6. 

Although in typical use, IDCOM and WRTID would be used in connection with 
the peak finding routines PEAK etc. (p. 97) using values of ES obtained from it, 
a small main program IDTST (p. 98) exists which takes data from cards and provides 
a summary of isotopes possibly present. 
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PF«lt HO, 1 



"Pftk NO. 2 



PFAK NO. 10 



"Jia?" KAIN PEAK 

0239 ^'*^^ fE*K 

X-RAV ^KC^' Fb 



PCAJC _NtU 17_-. tN£F£y_. 21i6._C_j((tV 

3 CL38 ffiit^ FEJK 



PE«_K NO. __La __ _tNERG* 22^,2, Kt V 



SIhGLE ESCAPt PEAK 







20 It-2J3 

__Ar_«iy...f«c»' FB 

tNEftcr ytb.O KtW 

9 EU15Z 
10 eUl52f' 
12 G0159 
15 HFiai 

ifl pn97^ 

25 aAUl 
1 NA24 


I^AIN PEAK 


PEAK 
PEAK 

"peak 

PEAK 


NO. 20 
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ENEfiCY 2667 


iC_K_Ey_ 
re KEV~ 


SfALL PEAK 

sue PEAK 
SMIL PEAK 








PF«X NO. 


3 


KAIN PEAK 




NO."^ 
NCV 
NO."" 




5 

5 


MN56 
MN56 

6RGY^"y753 
NA24 






fllt^ PEAK 

^AI^ PEAK 

fUS PEAK 
filh PEAK 
KAIN PEAK 
CCuBtE ESCAPE PEAK 

P-AU PEAK 






?1 
3 


EK 
1 


fAIN PEAK 










4 


ENERCV 511. C KEV 

8 SEe3 
17 CS190f 
19 TL202 
22 RUIQ6 
16 W187 






ENERGY 2962 
5 HN56 

e'nerg'y"" "aecT 

3 CL38 


.0 KEV 
C KEV 


S^ALL PEAK 










fAIK PEAK 
fAU PEAK 
(-Alh PEAK 
ifUL PEAK 
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PflK NO. 
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23 SR92 
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CCLeLE ESCAPE PEAK 
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6 
7 
8 
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10 
U 
12 
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10 
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27 

28 

50 

150 

151 
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CL38 

K«2 

MN56 

CC60 
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SE33 

EL152 

EU152H 

EU1S5 

GD1S9 


-- 


i- 




1368.4 
1778. S 
1642.0 
1524.7 
E46,9 
1332.4 


2753.9 

0.0 
2166.8 

0.0 
1810.7 
1173.1 


0,0 
0.0 

c.c 

312.9 

2112. a 

CO 


J.r 

CO 

c.c 

1923.0 
L.C 

O.o 


C.C 
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717.\ 
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■J.C 
22b. d 


0.0 
79S.7 






6 C060 


filK PEAK 


1315.0 




8 


ENERGY 1332.4 KEV 
6 CC6C 


CAIK PEAK 


58.4 


PEAK NO. 


13 
14 
15 
16 
17 
18 


158 
U6 
176 
182 

lee 

196 


TB161 

Tf-UO 

HFiei 

H1B7 

CS190K 

Pri97f' 








74.9 
84.4 
482.2 
685.7 
616.4 
346.3 


49.1 
O.J 
133.1 
479.3 
186.7 
0.0 


CO 

345.7 

72.3 

361.2 

C.C 


art.l 

CO 

4 70.0 

134.3 

50 3. C 

0.0 


0,0 

0.0 

615.0 

616.1 




9 


7 C060H 
ENERGY 1368.0 KEV 


CAU PEAK 


5i0.C 

CC 


prw Nc. 


19 
20 
21 
22 
23 
24 


2C7 
209 
211 
214 

233 


-TL2C2 
TH233 
U239 
RU1C6 
SR92 
1134 








440.2 
66.6 
74.7 

511. 9 
1363. 3 

846.7 


0.0 
459.2 

0.6 

622.3 
241.6 
803.7 


521.5 

169.3 

0.0 

1050.5 
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0.0 

J.O 






1 NA24 


CAIK PEAK 


1142.0 
155.1 




25 


241 


8Ai41 








189.8 


276.9 


304.4 


344 .J 


457.9 



ENERGY 1524.7 KEV 



K42 r'AIK PEAK 
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PEAK PiC. 12 
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5 MN56 fAIN PEA 



PEAK NO. 16 



ENERGY 2112 .0 KEV ' 



Output from IDCOM Test Program IDTST. 
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Card input is evident from the FORTRAJI listing (p. 98). The first card con- 
tains the "cooling" time in minutes. Subsequent cards contain one gamma energy 
each (in KeV) . Up to 100 energies may be entered. The last card in this series 
should be blank. Sample output is shown on page 54. Each peak is first listed 
separately with its attributes. Thus peak 3 has seven possibly contributing 
nuclides. The larger number of attributes at low energies is in part due to the 
assvmiption of a constant peak FWHM of 26E = 4KeV. A form used in later runs was 

&E = 0.7 + .OOIE (KeV) 

IDTST provides a summary table in the form of a list of elements found, with 
the number of main peaks in the library compared with those found in the actual 
spectrum. The print out also includes the main gamma energies in the library. 
This table is helpful in determining which elements are really present. Generally 
speaking if not all main peaks are present for a given isotope its presence may 
be doubted. 

F. Non-linear least squares routine . 

Let the set (Y,,X.)j=l, NbeN observations on Y and X which are presumed 
to be related by Y. = f(X., S'^..,S> ) where the 3.'s (i = 1, M) are M adiustable 
parameters. The subroutine NLLS determines a parameter set g., starting from 
initial guesses for 3 (say 3.), which minimizes the residual 






J*/ 

with respect to each 6. and thus produces a "best" fit according to this least 
squares criterion. That is, a set 3. is found such that -r-r— = o for all i. W. is 
the weight attached to the j observation. If there is more than one parameter 
set which minimizes R, the solution obtained by the program will depend upon the 
starting point, i.e. the guesses for 3.. Moreover the functions must be bounded 
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in 3-, In order to insure convergence. Despite such limitations the routine is 
of great usefulness and power. 

Its application to the analysis of Ge(Ll) spectra is two-fold. 

1) Energy Calibration 

2) Resolution of Doublets 

The algorithm by which the minimization is accomplished is described briefly in 
Reference 24. Since a solution cannot be obtained in closed form generally, the 
equation is "linearized" by calculating small increments , 6 . , to initial values 
of 3?, which will tend to minimize R, that is new values of 3.j 

are calculated for which A^ 



f(X., &]...)Z f(X 
3 -*- 



,,o...).y4i i 



such that N j!|\. 



X"i^-Ko-Z|-^;^'= 



^\ ,... 



A solution of these M linear equations for each 6. allows one to correct 3^ in the 
direction of minimum R. The new values of 3., 3., constitute new guesses for 
3. which themselves may be corrected by calculation of a new & . This iterative 
procedure may be repeated until a desired accuracy in calculating each 3^ is achieved. 
That is, until . 

4 i 



for all i where e is small. 



f^\ ^ 



-57- 



Equation can can be seen to be a matrix equation. Defining F. = (Y, - 



■^ - f(X^, B^...) 



taking the partial with respect to 6 , gives, 



N 




M 



W, ( F, - 1 






•= o 



j = i 



which gives 



X"- 



^ 






y ^uali-^is 



'ap^^^K 



^4K^^ 



}-' 



>f5.« 



i'^' 



where A is a symmetric matrix with elements given by 






^-< 



symbolically the solution for 6 . is then 



k = C' Yk 



where A is the inverse matrix associated with A . The standard error associated 
li^- iK 

with each parameter g . is given by 



^.- 



J 



A/^M 
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The arguments of the FORTRAN listing for NLLS are adequately described by the 
comment cards (p. 99). Matrix inversion is accomplished by the subroutine INV 
called from the UM function library. The functional relationship f(X,, g^...) 
is supplied by the user as the subroutine ARC. The fortran listing of ARG on p. 99 
is for a power series, 

f(x , e^.-'S^) = ^1 + 3o X. + 3, X? 

jlm 12j 3j 

Where U is the value of X, for a specified J. ARG must be called for each data 

point. The dummy, P(l), contains the number of parameters = the polynomial order 

+1. V is the value of f(X.,g) returned. The current values of B are stored as 

G(I), 1=1, M, and the derivatives, t^ , are returned by the subroutine as D(I) 1= 

J 
1, M. Thus a linear least squares fit would be obtained by using NLLS-ARG with 

P(l) = 2. 

Since the correspondence between the energy of a detected gamma ray and the 
position of its peak in the analyzer spectrum is very nearly linear for the 
Nuclear Data 2200 unit, the addition of a small quadratic "correction" term suffices 
for most analytic work undertaken on this instrument. Thus a value of P(l) = 3 
is customarily used but the capability exists of making higher-order fits for 
gamma energy determinations. Experience has shown that with reasonable initial 
guesses (within 20% of correct values) the quadratic fit using up to perhaps 10 
lines, converges in 2 or 3 iterations with e = 0.1%. 

For purposes of energy calibration less powerful techniques of course, may 
suffice but, in addition, NLLS provides a means of analyzing doublets. Given 
the functional form, Gaussian or skewed Gaussian, of each of two peaks which 
are superposed NLLS can estimate their contributions under a broader range of 
conditions than may be done well by eye. Some experience has been gained in 
application of this program but its need in analysing typically encountered samples 
is as yet not great. 
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4. Composite Program 

The program SPAN (FORTRAN p.XO) is designed to obtain weights of elements 
present in samples deposited on various substrates and irradiated according to 
the analytical scheme discussed in Chapter 2. For each of the four combinations 
of irradiation and counting modes, the program is supplied with a list of gamma 
energies for which it estimates the position of corresponding peaks in the analyzer 
spectrum at hand and calculates net peak areas. The general features of the pro- 
gram sequence can be seen from the abbreviated flow diagram (Fig. 13.) A sample 
input card set is shown on p. 61. Data on these cards were developed and actually 
used in analysis of air particulate samples. Up to 50 gamma energies may be con- 
tained in each of the four blocks of isotope data accepted by the program. Two 
cards are required for each gamma energy entered. In the program the four run 
modes are indexed as follows : 

ITY Irradiation Counting Live Time 



1 


5 minutes 


400 sec 


2 


ff 


1000 


3 


5 hours 


2000 


4 


II 


4000 



For the short irradiations (ITY=1,2) the isotope cards have the following format; 
NS E SM FKl DFl FK2 DF2 



CARD 



2A4 2X F10.5 F10.5 F10.5 F10.5 F10.5 F10.5 
POL DP SOR DS POR DR 



F10.5 F10.5 F10.5 FlO.5 FlO.5 F10.5 
NS = the isotope symbol (e.g. CA49) (ITY=1) 
■£ = gamma energy (KEV) (e.g. 3083.0 KEV) 
SM = weight of element (i.e. isotopic mixture not the individual radionuclide) in 
the standard. (Arbitrary units e.g. 99.8 pg) 

FKl = the factor converting net peak areas (counts) to element weight in units 
of SM based on some standard flux measure. (Say unit: 0.211 yg/count) 
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DFl = the standard error in FKl. (.0275 yg/count) 

FK2, DF2 are similarly defined for a second counting arrangement differing 
from that used above. (e.g. FK2 determined for a larger distance between the 
sample and the detector) . 

It will be seen that as an analytical technique is developed FKl and DFl 
are not initially known. If the program is used with those values omitted, the 
output will consist of unconverted net peak areas. These areas can be used to 
obtain the FK's which can in turn be entered in the appropriate spaces on the 
card. The computer could have been programmed to make this calculation itself 
except that (1) once determined, FK's should not have to be re-computed, and (2) 
the estimate of DF should involve not only the Poisson-statistical uncertainty 
of individual estimates but also the uncertainties in reproducing the value of 
FK experimentally. Thus FK and DF should result from a set of measurements for 
which are most conveniently analyzed once by "hand" using the computer to provide 
net peak areas. 

The second of the two isotope cards contains substrate data. 
POL = the amount of element in units of SM present in the substrate on which the 
unknown is deposited per standard substrate area or volume. (i.e. yg Ca on a 
full-circle polyethylene collector sheet.) 
DP = the uncertainty in the value POL. 

The pairs SOR, DS and FOR, DR are similarly defined, so that the program as 
presently written allows substrate or blank subtraction to bemade for three 
different materials. If no data is entered on this card, the program will operate 
properly assuming zero for blank substraction. It will be seen that given 
FKl or FK2, the activity of the blank material can be found initially by running 
the program for blanks as samples with zeros for the blank data. The weighted 
averages and standard deviations resulting from measurement of a number of duplicate 
blanks and analysis by computer, enables these zeros to be replaced. 
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The order of isotopes within a block is generally unimportant but in its 

9/: 0-7 
present version the program corrects the total counts from Mg (n, )Mg 

27 27 
activity (E=1014 KeV) for the interference arising from the Al (n,p) Mg reaction. 

In order to make the correction properly in the present program, the Al isotope 

cards must appear third in block with ITY=1 and the Mg isotope cards must appear 

sixth on the block with ITY=2. Moreover the analysis of the ITY=1 spectrum must 

preceed analysis of the ITY=2 spectrum. 

For the long irradiations (ITy=3,4) the isotope cards contain the following 

data: 

„. T T NS E SM ASH DA TAU FMY 
Sxngle 1. 



^ , 2A4 2X F10.5 F10.5 F10.5 F10.5 F10.5 F10.5 
Isotope 

^^^ 2. Same as card 2 for ITY=1,2 

NS, E, SM have the same meaning as before. (Say Zn 69M (ITY=3), 438.7 KeV, 
5.055 yg). 

ASH = the amount of element present in a standard-sized portion substrate 
containing the standard. (units of SM 0.1 yg) 

DA = the uncertainty in ASH (0.025 y) 

TAU = the isotope half life in hours. (13.8 hr) 

There are two choices for standards with ITY=3,4. Either SM or SM*FMY. 
The choice is determined by a parameter in the control cards described below. 

The ordering of isotope cards is again unimportant except: 

1) For ITY=3 the first card must be NA 1368.4 KeV and the second CU 511.0. 

The program corrects the observed 511 KeV counts arising mainly for the Cu g+ 

24 
annihilation for a small contribution arising from the energetic Na 2753.6 KeV 

gamma ray. Experimentally: 

Net CTS Due to Cu = total CTS - 0.093*COUNTS in 1368 KeV Na peak 

2) For ITY=4 the first card must be SE75 264.6 KeV while the third must be 
Hg203 279.1. Since the 279.1 line is the only one available for Hg 203, it must be 
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used if amounts of this element is to be measured but the, 279.6 KeV Se75 line is 
virtually coincident with it. This interference can be removed by subtracting 
the counts inferred to be due to Se75 based on the ratio of counts in the 264.6 KeV 
and 279.6 KeV peaks for pure Se. 

Experimentally Se(279.6)/Se(264. 6) = 0.385 so that counts due to Hg = gross 
counts - 0.385*CTS from Se 264.6. 

These interference corrections which presently necessitate ordering of the 
data and spectrum processing in the case of Mg-Al involve minor modifications of 
the program and can be removed if calculation of weights of these elements (Mg, 
Hg and Cu) are not of interest. 

The end of each of the four isotope sets is defined by two blank cards. 

The isotope data cards are followed bjr four corresponding sets of gamma energies 
(ITY=l-4) for use with an internal energy calibration routine. The list consists 
of energies corresponding to recurrent, prominent, isolated peaks. Energies entered 
may not necessarily be those in the isotope data set. While gross changes in the 
calibration must be treated by supplying information externally from cards it 
has proved useful to compensate for small alterations in calibration (shifts of 
up to several kilovolts) by means of information in the spectrum itself. Expected 
peak positions calculated using the above lists locate an interval over which a 
search is made for the absolute maximum. If the maximum occurs at the ends of 
the interval, it is considered out of range for inclusion in the final list of 
peaks for calibration. If it is within the range its statistical significance 
is checked. If it qualifies statistically, the centrold is located using a parabolic 
interpolation routine XM. This routine has the effect of adding the correct decimal 
point to the channel with maximum counts. The subset of energies and identified 
centroids satisfying both position and statistical criteria is then used to update 
the calibration prior to spectrum analysis. If the number of energies left in the 
subset is less than the number of coefficients in the fitting polynomial (order+1) , 
no update is made. 
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Once established, the above set (isotope data and internal energy calibra- 
tion data) is included routinely with every tape analyzed. Cards specific for 
processing a given tape follow this set. 

1. Label Card LABI, LAB2, (2A4) label of the tape being treated (e.g. MT27) 

„ ^ ^ n NFW NDEL NBASE NCUT LWD SFLUX 

2. Parameter Card -^4 J4 ^4 J4 14 110.5 

For each gamma energy an interval + NFW channels centered on the calculated peak 
position, ICH, is searched for the absolute maximum channel, IMX. The net peak 
area is calculated over an interval + NBASE channels centered on the IMX . As in 
the case of the energy calibration, the maximum must be within the interval ICH+NFW 
to be considered correctly located. If the maximum occurs at the ends of the interval 
a warning is generated. 

A measure of the statistical significance of each net sum over the interval 
IMX+NBASE, is computed and printed out. The statistical criterion is defined as 



CRIT=((N-B)// B ) where N = total counts over the interval and B = the corres- 
ponding baseline sum. Thus CRIT is the square of the ratio of net peak counts 
to the standard deviation in baseline counts in the interval. This expression 

for CRIT can be related to the decision, detection and quantitative determination 

(25) 
limits L , L and L defined by Currie in a careful discussion of signal de- 

tection thresholds. L is the value of that net signal which can just be detected 

in a single observation. L is the value of that true limiting mean signal which 

can be reliably detected in a given analytical procedure and L is that level 

at which an individual signal is likely to be close to the limiting mean (somewhat 

arbitrarily taken as being within 10% with a 95% probability). "Working" expressions 

(25") 
for the three limits are given for both Gaussian and Poisson signal distribu- 
tions. For the latter case 



L^ = 2.33 /T^ 



L^ = 2.71 + 4.65 /yT 



L^ = 50 (1 + / 1 + M„/12.5) 
Q B 
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for "paired observations" of signal and background where y = the limiting mean 
value of the background = (T^. For sufficiently large y , 

4.65 /lA 50 



fey ■".&>*■"&> 



/ 12.5 



The approximations are correct to within 10% for yJ5-70, 34, 1250 counts respectively. 

In the present case the net "signal" is N-B while the limiting mean background yg=B 

and ^=/~B~. Thus 

L 
CRIT. = (-^ )^^ 6 
^ ^B 

^D 2 
CRIT^ = (-—) Ai 20 

CRIT^ = (-^)^'«200 

^ 'b 

From examples of isolated peaks having different values of GRIT, shown in fig. 14, 
it can be seen that there is qualitative agreement with the above results. For 
GRIT 45 there is essentially no indication of a peak. The values of a are also 
iheasured and indicated, where 

Maximum count in the scan interval - B 



where B = (B/number of channels in the sum) = average baseline count per channel 

in this region. Here a^l indicating that the maximum peak count is comparable 

to the background fluctuations. For GRIT -^ 10 there is a detectable peak but 

little of quantitative use. The peak maximum is typically two standard deviations 

above the mean background (a*'2). For 20 " GRIT "^ 200 detection is reliable but 

fractional uncertainties in the peak area determination exceed 10%. These results 

are summarized below. The values entered apply with y > 500 counts. 

B 
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Qualltative meaning of the statistical measure GRIT 



GRIT 


a 


Gomment 


-^5 


^1 


No detectable peak 


5-20 


^2 


Detected with marginal 
reliability 


20-200 


3-10 


Reliably detected but poor 
statistics 


>200 


^10 


Statistics better than 10% 



NDEL defines the interval for the internal energy calibration and for the 
titanium (320 KeV peak) flux monitor peak search. With iKeV/ channel calibration 
NDEL is typically taken as 5 channels. 

NCUT = the number of channels out of the actual spectrimi which will be treated 
by the program. Entering the minimum value consistent with retrieval of desired 

information will decrease running time. Gamma energy entered in the isotope data 

41 
set with peaks in channels beyond NCUT will be deleted. (e.g. the 3083 KeV Ga 

peak should be analyzed. The baseline should be calculated correctly 2£+l= 

channels beyond this point so NCUT is taken as 3200). 

LWD = the moving average half-width for the baseline calculation (BASE) 

(e.g. LWD=£=30) 

320 
SFLUX = the measure of reactor flux, in Ti KeV gamma counts, to which 

the values FKl, DFl, FK2 , DF2 are referred, (e.g. 6000 counts. This number was 

obtained from computer analysis with NDEL=5). 

3. The main control cards containing coded instructions for processing individual 

spectra appear next in the sequence. 



iTG Inrec nbg ity nbe nsk nrw nbks ngal 



14 I 812 

ITG - tagword sought 

NREC = 1,2,4 coded channel number 1=1024, 2=2048, etc. 
NSK >0: skip NSK records 



FTR GONV TIME 



3F10.3 
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NBKS >0: backspace NBKS records 

NRW >0 : rewind the tape 

ITG should not be entered for these operations. One instruction per card 
should be entered. At the end of execution of the Instruction, the program expects 
another main control card. 

NCAL <0 : read energy calibration from cards in the following format : 



EG 



2F 



POS 



10.5 



EG = actual gamma energy 

POS = observed peak position 

The order of the fit is - NCAL. Up to 100 such pairs may be entered. The 
end of this data set is indicated by a blank card. At the end of the calibration 
calculation, the program expects to read another main control card, so ITG should 
not be entered when NCAL<0. 

NCAL>0: after reading in the spectrum with tagword ITG an internal energy 
calibration of order NCAL is made before analyzing. 

No-operation values for the above quantities are thus zero. A totally blank 
card indicates the end of all data processing. 

ITY = run type index defined earlier. 

NBG = code describing the type of processing required for the spec,trum. 

NGB 
12 3 4 5 



ITY 1 
2 
3 
4 



Sample Spectrum 
Type of substrate 

is 
POL SOR POR 



Flux Monitor 
Spectrum 



Standard ■> 

Spectrum 

No Mulitplier 



Reset Flux 
correction 



Standard * 

Spectrum 

With Multiplier 



*Blank is ASH 
NBE =0 use FKl, DFl with ITY=1,2 not used with ITY=3,4 
= 1 use FK2, DF2 
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FTR = fraction of substrate (POL, SOR, FOR or ASH for the standard ) actually 
used. 

CONV ■■ divides the calculated weight of sample to convert to desired units 
(e.g. CONV = no of cubic meters of air sampled in which case output 
is in yg/m ) 
If ITY=3,4 and NBG=5, that is for a standard spectrum with modified 
standard values CONV = the multiplier of SM for those values of FMY in the 
isotope data set which are left blank. 

For values of ITY=1,2 NBG=4,5 
and ITY=3,4 NGB=4 

CONV is not used. 

TIME = time in the format hours -minutes for the onset of counting either the 
standard or the sample. This value is needed only for ITY=3,4. 

The program calculates the correct fraction in hours from this value. Since 
the standard and the sample cannot be counted simultaneously a correction must 
be made for the decay of one relative to the other. Only the relative values of 
time are needed in the calculation. 

Sample output. Computer output listed on pages 73 to 75 results from the 
main control data card sequence listed on p. 61. The first card with NCAL=-1 
produces a linear external energy calibration prior to reading any spectra using 
the energy cardswhich follow. The associated output gives the number of iterations 
required to produce a fit with e=0.0011, the linear coefficients with the inter- 
cept = -0.207 channels in this case and slope = 1.0037 channels/KeV. The five 
output columns are 

E C C AC CRIT 
o c 

E = isotope gamma energy 

C = observed peak position 
o 

C = peak position-calculated from the best fit 

AC = C -C 
o c 
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In the case of the external calibration CRIT is undefined and taken to be zero. 

At the end of the external calibration a card with NSK=8 is read and 8 records 
are skipped. The spectrum with tagword 9 is sought and is processed as a flux 
monitor spectrum. The associated output for this analysis includes a printout 
of the main control card image, the current energy calibration parameters and 

ICH = the expected peak location (319) 

IMX = the observed maximum in the scan interval (320) 

"FLUX" = the net peak area of the 320 KeV Ti gamma (6970.4) with the 
associated Poission-statistical uncertainly (89.5) 

NBE = index of the set of conversion factors to be updated (in this case FKl) 

SFLUX = the net peak area on which the conversion factors entered on the 
isotope data cards are based. 

FL = the ratio SFLUX/FLUX 

The flux monitor spectrijm for the short runs (400 and 1000 sec.) must be analyzed 
first before treating the sample spectra. In typical measurements the flux spectrum 
is generated after the 400 second count and before the 1000 second count. Thus 
in order to analyze the 400 second spectrum the computer is instructed via the 
next card with NBKS=2 to backspace two records. Spectrum with tagword 8 is then 
sought and treated as a 400 second spectrum.. (ITY=1) 

Here NBG=2 indicating substrate type tvro (SOR) was used but only a fraction 

equal to 0.164 of the amount entered in the isotope data cards is to be subtracted. 

3 
CONV = 0.0053 = .001*5.3 m of air were sampled by the substrate section and 

3 

the output units are to be expressed in ng/m while the conversion factors are 

in yg/count. (S is entered in mg/count so that the output units for this case 
IS yg/m ). 

Following the listing of the control data in the output for tagword 8 are the 
results of analysis for each isotope present in the isotope data set within the 
range of NCUT. This table has the following entries 
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IND, ICH, IMX, NS, E, CT, DCT, W, DW, PCT, FC, CRIT 
IND = the list index 
ICH = expected peak channel 

IMX = absolute maximum in the scan interval ICH+NFW 
NS = element sjraibol 
E = gamma energy 
CT = net peak area (raw data counts) 

DCT = Poission statistical uncertainty in CT 

3 
W = the corrected and converted net peak area (here in ng/m ) 

DW = the combined uncertainty in this value 

PCT = the percent of substrate in the total sample weight 

FC = the decay correction factor applied by the program. For ITY=1,2 FC=1.0 

CRIT = the statistical criterion defined earlier. 

It is assumed that the live time operating mode handles the dead time correction 
properly and that short-lived isotope dead time effects not correctable by live 
time operation are small. This is likely to be true if the overall dead time is 
kept low (e.g. <20%) . Also it is assumed that the short-lived counting is for 
the same live-time period (400 and 1000 sec) as for the measurements of corres- 
ponding FK's. If not the correction must be included in the value of CONV. 

W is given by 

W = (CT*FK-FTR-BL)/CONV 
where CT, FTR, and CONV are defined earlier. FK is either FKl or FK2 and BL is 
the activity of the appropriate blank, POL, POR, etc. The error in W is given by 



DW - 



(S!""'- ^ ©••"'' *(t J'""' 



where DEL = DP, OS, etc. 
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For the first Isotope, S37 (3102.4 KeV) , the peak is not found within the 7 
channel scan interval centered on channel 3112 but at the end channel, 3109. The 
net sum in the interval centered around 3109 is evaluated but the remaining cal- 
culation is deleted and a warning message printed. Inspection of the value CRIT 
suggests that the value is of doubtful statistical significance. This can be 
seen further from inspection of the original data which is presented below the 
isotope table. The contents of 23 channels centered on ICH are printed-out so 
that computer results may be readily checked by hand in those cases where any 
doubt about the analysis exists. The form of this table is 

1. Index line 0-22 

2. IND, IS, COUNTS (IS to IS+22) in 15 FORMAT 



Thus in the case of S37 with IND=1 IS=3101, the column labeled 11 corres- 
ponding to IS+11=3112=ICH, is the expected peak position. 

Upon completion of analysis of spectrum 8, 48 records are skipped and spectrum 
59 is sought. This is treated as a long-irradiation standard spectrum (ITY=3, 
NBG=4). Comparison of the card input with the control line print out (p. 74) in- 
dicates that that time entered as 9h.l5 min. has been correctly converted to 9.25 hrs. 
Also a new energy calibration update has been made which is not indicated in the 
write up,, hence the new linear calibration coefficients - 0.522 and 1.002. Output 
for the standard spectrum has a slightly different form: 

IND, ICH, IMX, NS, E, CT, DCT, FKS , DFS , PCT 
where all entries have been defined previously except 

FKS = calculated conversion factor (pg/count) 

DFS = Poisson statistical uncertainty in CT combined with the blank uncertainty DA. 

FKS = (SM + ASH*FTR)/CT 
and DFS = 
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and it is assumed that there is no uncertainty in the value of SM or FTR. 

It should be noted that while DFl is based on a series of measurements in- 
dicating the overall reproducibility of FKl, DFS is based on a single measurement. 
While it includes the established uncertainty in the substrate activity (DA) it 
is not strictly analogous to DFl. 

PCT = the percent of ASH in the total activity = 100. X ASH/(SM+ASH) 

With analysis of the standard spectrum completed, one spectrum is skipped 
and spectrum 61 is analyzed as a sample spectrum (NBG=2, ITY=3) (p. 75). At 
the end of this analysis the job terminates. 

Further information on all programs can be obtained from Dr. John Robbins, 
Great Lakes Research Division, I.S.T., North University Building, Room 1112, 
Telephone 764-2420. 
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Llst of Figures 

Fig. 1 Measurement of Linearity of the Gamma-ray spectrometer. 

Fig. 2 Determination of a net peak area. 

Fig. 3 Irradiation-counting schema for short-lived isotopes. 

Fig. 4 Gamma-ray spectrum obtained after 3 min. cooling time. 

Fig. 5 Gamma-ray spectrum obtained after 15 min. cooling time. 

Fig. 6 Irradiation-counting schema for long-lived isotopes. 

Fig. 7 Gamma-ray spectrum obtained after 15 hours cooling time. 

Fig. 8 Gamma-ray spectrum obtained after 15 days cooling time. 

Fig. 9 Flux gradient at irradiation site in irradiation bottle. 

Fig. 10 Ge Li Spectra General data reduction scheme. 

Fig. 11 Dependence of the calculated baseline on the moving average half-width, e. 

Fig. 12 Gamma-ray spectrum showing the calculated baseline and computer-located 
peaks. 

Fig. 13 An abbreviated flow diagram for SPAN. 

Fig. 14 Sample spectrum peaks versus the statistical criterion, GRIT. 
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Irradiation and Counting Scheme for Short-lived Isotopes. 
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Irradiation and Counting Scheme for Long-lived Isotopes, 
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Fig. 10. A General Data Reduction Schematic Diagram 
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An Abbreviated Flow Diagram for SPAN 
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The Main Tape Inventory Program INVEN : 



th IV r. CQ^PILE" 



M* IN 



C3-Zfl-TT . 



p><;e nnoMN iv g crwpiLCR 



■IAIN 



0^-'?9-TP 



C "tESSY T»PE «E«0. INVENTrRY »NC SUPER DUPE' 
C T TRACK BCO RE»D AND WRITE 
C R6*0 FROM PSD^ 2 WRITE UN PSDN 4 
INTeGF««2 LEN 
PIWEKSION ISIft^f"! . 



CLEAN WRITE 



CALL RFAn(ISTLENtl'iLN(4,etin,eni,C8l,Cait.CI!l*C0t) - 
t*Ftt»lFN 

HEN»(»EN/6I-? 

WRITE (6.1001 J,ISt3l,LEN,<4E^.J 
1^0 FQBMATIl^^K.'FILF ^0. ' . M.' TGMO *.tA*,' tWGTH- 



Cir'FftSttK LA(2')t 
0!"»EKSI1N IBtZi 
DATA H/Z^COOOOOO/ 
READ l^.in^l LB( 11,L6<?) 
104 F0RNAK2A4I 

READ I ^.Tm I ILAtJNJ - UgCt 



fl'KLS •.!*#• CU** RCS 
GO TC U 

525 FORMAT^f.* END OF FHF • 
GO TO a 



• 141 



700 F0RKAT(?0A4J 

REAO 15.411) N.KSK 
400 F0R«AT(2r4t 

IFINSK)6t6<T 
7 cnNTIWJF 

C ALL SKIPfl.NS K,2I 



526 FnPMAT(5lt 
a CGNTINIIE 

517 FORNATIIHI I 

CALL HEWIMO (*» 
STOP 



THERE IS AN ERROR ON THE »4EWLY WRITTEN TAPE 



6 CONTINUE 

RSF--0 
KNTF-fl 



RtRB-n 
KPEt-<^ 
KCUM"fl 

URITE(6.60II 
601 FORMAIIIHI ■27K. 
1 CONT IMIE 



CONTENTS OF TAPE BEING BEAD *.////» 



101 



RCUP-rCU*"*! 

CALL PFAOnS,tFN.M,LN.?.C10,611,Cl?.ei3,£14,SI51 

••EH-LEN 

KF«KF*1 

WRI TE 16,1011 KF . I SQl.H N.MEN.KCU'* „ _, 

TnHMATI34X,*FLN '.U.' TGW •,IA4,» LEN '.I*.' CHNLS ' 



CMBO ' 
KNTF«fl 
IF(KeN-ir?4l31,21,?2 , ._ _ 

22 lF(MEN-2n4S)31.2l.21 

23 IFtHFN-40q6m>21,31 , 

2V WEC*K«tC«i 

WRITE (6,5001 ((OEC ._ - 

500 FntlMATnH*t«5»»' CRN ', I3l 

CALL KRITEIIS.LEN.n.LN.41 

KWT«'< „„ 

1FtMVLMSI3t)-NIU3,3 
WL IS FUNCTIO** TO CONVERT A TO I t«OOE USING BVTE-WfSEJHIFTIMG 
1 CALL REWTN0(4t 

WRITE (6,524) LB(1I,LP(2> __ — 

524 F0RP»ArllHr,*4ir, 'CONTENTS OF NFWIV WRITTEN TAPE S2A4,//I 

W R1TE(6.B001 ILAtJ),J-l,201 



2 CONT!KUF 
GO TC 1 
10 KEOF-KEOF*! 
KNTF*KNTF#1 
WRITf 16,2001 KECF 
2nn FORHATIZOX.'EOF ', 



14) 



IFIK^TF-nll.I, 
5 WRITE(6,513I 
503 F0PMAT(1H*,30N,*SK|P*I 
CALL S«IP(0,1,2) 
KNTF-fl 

KNT-n 

GO TO I 

11 WC«fl 

GO TO 16 

12 HC-12 

GO TP 16 

__ 1 3 HC-16 

GO TC 16 

14 nZ'itt 

GO TO 16 

15 MC-24 
WRnE(6,l'\'13i 

10^3 F QRMATI20X,' TAPE MACHINE FAILURE COOF 24 ■ 
STOP 

16 KFPR-KERR*! 
WRITE(6,360I KERR.MC 

310 F0RMAT(20X,'ERR ',14,' CD ',14) 
KNT-KNT^l 

IFIKHT-2tl,4.4 

4 CONTINUE 

CALL 5KIP(0,l,2) 
KNTF-O 



SKIP'! 

GO TO 1 

31 KSF-KSF*! 
IFIKSF-411,32,32 

32 WRITE(6,5J1) 

531 FORMATIJOXf'quEFR FIiE_SKIP'J_ 
CALL SKIPi0,l,2i 
KSF-0 



^L MEMORY REQUIREMENTS pOBCfrA. BYTES, 
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The RDLST and BASE Prograips 



_!tAU-_ 



,.,** yt s cf yuEit 



a«E 



C "*1N KE»C List PRCCBtt- 

CALL BFWthDt*! 
1 cnWTIHUE 



SURHOUTINC MSea.lST»INC.I»,8*l 

cni)>\Vt WECiSIOM •CM 

OOUBLE PSRCISIOW B$l*0»» . 

DnuUtt MECISIOM WUOMI 



ilLAl)(»*tf-)| NIOihPtCtKSKf IS(I£.hKS 

CALL SKtPtn,NSK,4) 
NP.KREC-t 

2 STC P 



3 CONHr-UE 

a cnhi IMUF 

CALL PCTAPEISP.n.r 

irc»Kvtini 

ir ILTC-MTG)fl*4<ZM 
21* WftlTEI6t215t LTCt HTC 



215 FCttNATI////.in]( 

srnp 

IFtMMSIllltlU.U^ 
in COHTIKUE 

C*LC WBTfTI IS, le. i3»n» i"i 



TACkfUO "E«C»»fl5,* I»Ct*DBC SOJCM- 



ACh-?»L»l 

IE»!fcD-L 
ISF-IS*l 

KK-IS*1 

DO 1 J-!5T,!WC 
I SSI J)*I"tJI 
on 2 J-l.!** 
AJ-IO-J 
FD-AJ*0.2 
W( IS 1-1.0 
DO ?S M»|ST,ISF 



"'GO TO \ 
112 CONTIMUE 

CALL eASk(r<&S>iS(IE<IPtBSI 
CALL MRIASIIS(IEtI9tl1.|RtSS) 

Gfl !0 I 

ESQ _. 



"25 MllSI-tlllSI^Stfl 

DO J *-lt»t,IE_ 

ll|KI«WHC-il*PlSlil*lt-ESIK-t-lt 
3 CnnTIHUE 
00 S K«i$,IE 
IFIMIKIIT.a.t 



8 CONTINUE 

MIKI-WIKt/ACH 

A-SSIKt 

q>W(Kt 

If tH-l»..FC*SWTIBn>,»,» 



l*0»Y HtCblPfHtKTS ' >ClHjnrUS 



5 CONTINUE 

7 CONTINUE 
BFTUHN 
END 



y.G CrWPjLEP . .__ _ ^JDT-S Pf 



!^-:7->_ 



SUHOnuTIwe oOTAPE(»l''»n.IT,IP) 
nH*E»i^tn"' jaC*-"!?! 
.2 r'>KTIMJF 

■FAm*,ni.f.*<']»it ii.tz.PtiT.iitiji.j-uNPi 

11-^ f<ll>'«ATI3A*.?;>l l? » , l 6>> 

BFTUP« 
I taftlTri*..?''"!! 

3nn FnPW«T(t EvU-fF-FtlF M 
GO TT 7 

EWn 



>N TEOWINATEO 



- Jt IV G CQOilCB 



MKTKS 



0>-«-TO 



SMtOUTI NE WRTBS f IS t IE* II* IT , !■*•$!_ 

OOUBtEmECISION »S(4(>4S> 

DIMENSION Iti4n»<n 

DIHEhSinw ICUTIIOI 

WPITEl&tlOl) . 



*| I* 6 CQWILEP 



n?7f PVtFS 



iSKL- 



I'fl FOKIiATlim » 

HttTE t«.941t 13 
" ■(•W* FOKKATI JO««*TACW)«D 
NW>tIE-IS*I)/10 
00 3 J-1. Nh 

js-if ij-n*i$ 



IT 



" ll»f Tim» •.I6t' SEC. 



Function "vliii 

IKfEGER 5MFT» 
L1-SHFTI((I.241 
Ll'LANOIllftSt 
L2«5HFTR||,UI 



WtITFI&«20f) 
:>nn FOftNATI/l 

imiTE<6t3OTI JS,ll»ltt.i-JS*Jf**JF 



L2-IANDIL2.19I 

L3-SHFrn( t.ai 

L3-IAN0I13*19I 

L4>LAN0I ltt9l 

)«Vl-l**n"(L3*l''»tL2*lC«llH 

PE TUPN 



e'nd 



K-L-JS*! 

IOUT(Kt>SSIL> 

iri C0NT»i*yE 

MtlTEfA.A'^TI linuTKI *K-t(m 
(.10 F0RNATI8K. inilX.tPII 

3 CPKTIKUE 



I mtnr^ PEOUIREHINTS OtIlFF BYTES 



I NEMORV «;EOUIREN£NTS OIOSE* BYTES 



M IV g CONQUER 



WRTW T . 



SU0RDUTINE HRTNTIIS,IE*l3,ITt im 
OlMENSION IRI4099) 
MITEI6,I234I 
1234 FDPNATIIMI I 

WRIT E 16.1(101 a, I T 

5M FORWATISOK.'TACWWD- ' H4.' LIVE 
NU*IIE-IS*U/10 
00 *3T J"t,NW 
JS»10«IJ-1>»IS 
JF-JS*9 
43T MBITE 1^. 40^1 J S.l lBltl .1 «JS j J F t .-Jf 

ttVt Foi»«»Tii»,r4.m.ifmJi,iei.4x,,i4i 

RETURN 
END 

I m'fmorv requirewents ooozcc bytes 

a — 



VlHHnilTlNF PKHf ( ICHTUFM,IR,HStRAT|n,SIIMtSO) 

DOUftLE PHECISION. SftS ^ . ._-„__„ 

onuRLF PHECISifW BS(*09S» 

OlKEMSlON IhUO^SI 

!<;»irH-NFw 



Sl.fc' 



=n, n 



U«IR( J-l 



IFim71il71.172 
17) KAT!o«-n,n 

GD..JQ-.113 

IT? f.lMtlHMP 

HATUl-SllM/T .,. . 
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A. Simple Main Plotting Program » 



FORTRAN IV G COMPILER 



f^A I N 



1 0- 2 4-6 9 



l5tQ<t.22 



PAGE 0( 



0001 



0002 
0003 



0004 
0005 



0006 
0007 



0008 
0009 



DOUBLE PRECISION BS(4096) 



DIMENSION SR<<i096),CH(4096) 
DIMENSION IPI4095) 



FTY=2.C 
FTX=.0O25 



XMX=66. 

CALL PLTXMX(XMX) 



XST=-ll.O 
999 CONTINUE 



0010 
0011 



0012 
0013 



READ(5,42C) NREC i ITG.NCUTt NSK.NCR 
CALL SKIP(0,NSK,») 



420 F0RMAT{5I4) 

IF(NRECI805,8C5,717 



0014 
_001^ 
0016" 
0017 



717 CONTINUE 
NR=NREC-1 



478 CONTINUE 

CALL RDTAPE(NR,I3,IT,IR) 



0018 
0019 



IFIMVLn3)-ITG)478,799,805 
805 STOP 



0020 
0021 



799 CONTINUE 

XST=XST+11.0 



0022 
0023 



XND=XST+10.0 

CALL PENUP1XST,0.5) 



0024 
0025 



CALL PENDN(XST,0.5» 
CALL PENDN(XN0,0.5) 



0026 
0027 



CALL PENUP(XST,0.5) 
CALL PENDMXST,9.0) 



0028 
0029 



DO I J=l, NCUT 
SR<J)=IR(J) 



0030 
0031 



IF ISR(J)) 2t2,3 
2 CONTINUE 



00 32 
0033 



SR(J)=0.5 
GO TO 987 



0034 
0035 



3 SR(JI=0.5tFTY*ALOGl0<SR(J)> 
987 CONTINUE 



0036 
0037 



U=J 
1 CH(JI=U*FTX+XST 



0038 
0039 



CALL PLlNEJCH(l),SR(l),NCUT,l,OtOil» 
IF(NCR)g98,999,998 



998 CONTINUE 

CALL 8ASE(30tl,NCUTtIR,BS) 



0040 
0041 



0042 
0043 



DO 51 J=l, NCUT 
T=BS(J) 



0044 
0045 



IF«T-1.0)52,52,53 
52 T=1.0 



0046 
0047 



53 CONTINUE 

51 SR(J)=0.5+FTY*AL0G10<TI 



0048 
! 0049 
r 0050 



CALL PLINE<CH(l),SR(l»,NCUTvl,0,Otl) 

GO TO 999 

END 
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The Peak Finding Routines 



'.^TP,lPK,<.lVi 






►■=1 ll-^O 

TALI UlllNni^l _ _ _ ____^ __^____ I2-0 

7 ccjrrjili . - - -- -- .-. - ^^^^ 

cm ^h if (ciSK,-!) ■ " " i^iij 

100 Ui-'MAT fni „__ T=:«r> 



1)0 I J=l, 
WKlTFIc? 



7 CCNil.'iUE 

""■ "" "" " ■ — "T-rrtU-n""" 

V CnvptLF s CIJhZ lO-.-'^i--.-) iMf-.'f, PAr.F n.,61 - - ''"--'•?i''.f.'M^l 

c ■ ■ " \z-f.0''WL8 tF'i'icRiTicVs ' ~ — .—— — -^_— — „ i) r.n to (i.i2.i3i.Lc 

c s^^cF iHE hHLiE t-'i..'v c rs iu,nin\fii c-c Toi 

c Ki< np= 10; .y< iTrTTir —" ij l"3=l 



OIMENSICn C(2j2).S(>C.?)" ^ 



C199l' l.u 11 II- I 



U GO lr> (14.1 ,l5l,Lf 



_Cl!iMl = l,_0 _ GO IV- 1 

fFff7;i(>r6.?'"^ " ■" ^ ■ "^ ' — " " ■■- !5 I3-I 

_2_ or 5 1^1, NZ MSI-1^_ 

" ""NT = i"*"^4.1 — • ' — " " "GC toi 

___M=iOJ-\T ____^ 1^ f^? te n Ife, 

Nf'«K..J*i<r ~~~ ~- ■■■ 16""16="I5" 

00 3_J = MS.f;F__ I5'I-I - 

■j"s-J-<< " ^ ~ — ' — — — —■ — —- iTTsT^TS 

IFIIIM, 



5U)')"" " " "" '""" ■ -• .•.o IFli2-U»l,el.hl 

DC i K 'JS.JF "■■ -.u I Ff I" jiTrjT7tT,5"2~- 

3 SIJ) = sniKlK) — — — ,^2 IFH5-[311.I,t3 

nc 5 J=^s.^f _ __ _ ^^ "■^■" 6 r"C dnTTnle — — 

5 ciji-siji " — - — u-^ Bs^s^ h(l^n 

-..-Cs-o , ... U=li 



• CS-CS*C( Jl*CtJ) 



~!^nrCI^TT^HJ^- 
NTST«19 



' ^ covPiLEH iut lL^nzh2 n•.?^.25 , p^gf Q.mi "" ■■" " """!!^!T: 



^UBKniTINE SFZWrS?w 



1DX-M7*"»1 """" To N^rS-TjTr^ 

IfN^ftUOS-IDX — NTSI' 20 

JS'ipJ(tl_ _____ _^___ ' TFHfl«S( N-'M)-NETiTT 

nil 1 j-is.iFN " "" : i\ N-n-12-i 

^SZM j|=c.C^ _ _ ^^ STSTiTI 

jF«j*i"rj)!" ■ """" "■" — — ~ — — 



ftc 1K-JS,J"F - — -^ !f irj?ITI.71,7; 



"^ L - ____ .^__^ '^ N^'i y 

[STkl " " — — "- J2 IF(M?-Wil,2«,2', 

" -'_!*f^_'.''.!<i'j . J* rrn^-in'fjs.fej.ft* 



"FTUPti 



N-r?-!!*) 



^ 2^ CDNTIMU! 

f lil-Jt.TIlA F K 7 1 K ; , 4 , f .-. ;^ ) SUMfl-tf'.C" ~" 

S=G>>'/<:".3"5'1 """ SUHT=P.{1 

IF((j7J1.i,; ____ "Po 2"6"~R=13,15 ~ 

■ M •?■* f'4 i " "■ " "" — SUHti = Stl>'P*S? VMKI 

T * w* w/ ?. . C AK^lf ~ 

"TSTTTT^TTTZ — — — ^ ^ ^^ SlJHT=Sti fn»AK « sZW 

s=5*Si,>:Tio.s*o.:)*';i:^Tii.o*T /(S"i))i __^_^____^^ iFiSLi"mpT,i."eT" 

. ccNtPiut "" (17 crtjTuur 

>'<Z-P^IJt';/GA''< _ _^„___ tPK-IPK»l """^ ~~~ 

'■f^TlJtN" ----.. .- -.- - -^ ... __ STRI tPKt-SLWT/SUI 



■Mnny BcuCi'iremi-nts rooAir i 
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Isotope Component Identification Programs : 



IV G CDHHIER HAIN 10-21-4S C9J5B.2* PAGE 0001 


IV C CflCPIlEB IDCCM 10-27-M CSiSe.27 PAGE OCCl 


c fith pncG rc AUK idccc foce c«T« on cjircs 


SIBRCL'TINE IOCO«Tl«,IS. I W.N.NI .ICX.KS.KPP. IP. INDX.NPKL ,ERI 


OlMhSION ESIlOO). NIICCI. Nll2C»100t . ECXUOaOOl ,NS1 lOO.^t 


OIHEhSION ES<100I.NIlCCIrkI(20,lCCI>IO]>l2C,100).KS(tOO*2I.IPU001 
OlFEhSIOK IhDXIlOOI.NPKLIlOOI.ERIlOC.SI 


READ 15.1001 TIME 

K-0 


OIHENSIQN KAHitt ECAPSI60) 
CALL DEWIN0I3I 



___^ BlAD(5,J0CiJ^(Kt 

100 F0ftNAf(2f 10,») 

lf lCSUII2.2,l 

2 IPk-k-l 



hLP"2*3 
D0_ l_J"lj IPK 
1 NtJl-d 
hPP»0 
DO 2 f0U*>-l*NlP 



CALL lOCCKITIME.ES.IPH.N.M.IDX.NS.NPP.IP.IhC^.hPKL.EPl 




READ 13.1001 I.N1.N2.TH. 


(ECAHIKI) K-1 


.91 


bll|TEU.201) 
201 FORMATtlKl t 




100 FaRNATti4,2A4.«F10.3l 
1F<TIHE-10.*TM1«!.2.2 






CO 3 J-t. IPK 
tilllTEI6.200) J, EStJI 




KK-0 






200 F0R«»TO0X./7/'->eAk *C. • Vl J.lOJt. 'EVlflGy ■»fT.l,' KEV 
CALL hRTlCU.N.NI.IDX.IP.MS.NPPI 


• /> 


E«tGJfl|KI 







CO * Jaf. APP 
_ Milt (fciiOOi ■J.lF(J I.I>S <J.U. RSIJ,2 1,l»IDXl J>,HPKLIJt, 
" 1 IERfj.K).k-t.5l 
300 F0BHAT11QX.M.2X.14,2X.2A4.2>.U.2X,H.512X.FT.111 



IFIEIT.7.6 
ffUr5l»-C13i;7i3_ 

iPR" 



3 00 4 J«I* 
If lABS<e-E StJU-EfieSI. I) ll5.5f» 











conriiiuE 








■ 


IFINMI.I 


DIMENSION Nl 1001. NII2C.100)t 10X12C.100I 
OlMEkSlON IPtlOOI. NSI100.2). 12(09} 


•SCS'II 






kSIIIF>.l)-lll 
NSINPPtZI'NZ 


DATA 12 /'M *.>K£'.*Ll>.'eE'.'8 '.'C '.•» <,* 
l.'MC».»Al»,'SI»,'P •.•S •••Cl'.'A '.'H •.•CA» 


,•¥ •••CR" 




NPKlUrF)-KCX 
M 57 «■!. S 


3.«SR»,»V •,'ZR».»Kl».«PC',»TC',*RU».'Rf''.'PD' 


•ACS'CO 


.•IN'.'SN' 


«7 
2 


ER^NPP.XI-ECAHIX) 
COHriMJf 


*.'SB«.«TE».'I •.•XC».«CS'.»flAS HA •.•€£', 'PR' 
S.'CO«,'TB».*DV».»KC',«ER'.«TMS'VB','LU',*ff» 


•MD',*PM 
'TA'.'K 

trn'T'Tp" 


.'SM'.'EU* 
.'RE*. 'OS' 
.'RA'.'AC'/ ~ 


— 


1 SO-0 

00 9 L"l» liPP 


NJ-NU) 


|N0X|LI*0 



"fMOxTi.Ji 

CO <> Ifl.fcFP 



i^U-IP(K}l9».99.«a 
»■ CpW TIWU E 



14 NIJ)-AUl*t 
. M11HIJI.JI»2 



-DElU4,14.1$ 



' "SJAtU'NSU.ll^ 



DOUBLE ESCAPf PEAK • ) 



203 MRITCU.MS) K, NStK. ll.NSIK.It 

M y FORMAT(2dXtI4.SXT;A*,5X,* SlRgLE ESCAPE I 

NSIKi ll.hi m.2l 



2«» MlirCI*.309l K. hSIK.l).NS(K.2l 
M» MRMATI20X,I*.5X,2A4.3X,« SUP PEAK 



CO TO 4 
; tO* COW T IMUE 



•X-RAV FRCP" '.M J 



60 ro 1 

9 hRlTI 14.4 171 



*PEAK hCI AI7RISUTECM 



■Q-ESiJIl-CELIIfc.lft.lT 

^X^i*S4i^Hil*l ___, 

NI(NIJI.JI-3 

_„ I CXtNUI.J»»l ' 

17 IFIK-4I1B.1B.I3 

It 00 20 KK.KS.5 

EE-ERU.KKI 

IFIE E)20. 20^21 

31 SUIM-E4EE 

If l ABSlSUIt-ESI JIHDEil22.22.20 

20 CORTlkUE 

60 TO 13 

22 NIJ)*M(J1«1 " 

____Nim(JI.JI*5 

lOXfSrjTnT"! 
13 CO NTIN O E 



READ (3.3001 (EGAMSIKI. 
300 FORMAT 160F10.3I 



_CUMI ICW EF(X1 



ICXlNULJIoI 
ji CO WT UUE 



IF(E5(J)-100.I29*29.31 
29 ET* 1 .132»ES(J1 



IF(£SlX|-iaO.I33.33.3C 
33 1F(AI)S(ES1JJI-ETI-1 .5I34.3A.32 
32 CONTINUE 

60 TO 30 



X-ESIJl 

I0X<MJ1.J>-29.033*.W»»X-.0030B5«X*X 

N(JJ)>N1JJI«1 

NI1N (JJ1,JJI*6 

IDXIN(JJI,JJ)-IOX(N(J|.J| 

)0 COWTIWUt 

31 CALL REWIND 131 

fiiiyRN . 
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Non-linear Least Squares Routines : 



SUBKaUTIKC M.LSIX.Y i >i.U .G.PiSE.VA B.E««.EPS.NlTiN D.NP.KI 

C NCh-LINEiK LE«ST SUU4R6S flTTlNG RUUTINE 

C LiXi PAIR S Xlll.Yll l CESCBIBEC B Y A CONTIN UO US tNO OIFFERENTI A BLE FUNCTIOW 

C F(X,G(JII xHERE GUI'S Ifi. ADJUSTABLE PARAMETERS. NLIS STARTS WITH GUESSES 

C FCR GIJ) <^C RETLRNS "BEST FIT " VALUES FOR THEN WHICH MINIMIZE THE 

C SUM OF THE SCUARES CF Il-E CIFFEHENCES BETWEEN Y AND F. 

C ^ 0« NC O F CAT APCINTS » (JI. Y(JI 

C W'STATISTICAL HEIGHT ASSOCIATED WITH YIJI. «1/STAN0AR0 0EVIATI0N**2 
_C GIJI«GUESS PA RAMETERS INITI ALLY, -BES T FIT PARAMS. UPO N RETURN FROM NILS 

C P»DUMHY PARAMETER CCMCUNICATIKC WITH ARG 

C SE(JI»STAKCARe ERRCR U PARAMETER GUI 

C VAR-VARIAKCE CF FIT 
JC ERR 'O UN L ESS PET OF UVERSE M ATRI X »0.0. TH EN ERR-I 

C NIT«NC OF ITERATICKS ALlOwffi" 

C K -AC TUAL KUM BEP PERFCRfEC 

C EPS-ALLOWEC FRACTICNAL ERRCR IN EVERY GUI TO TERMINATE ITERATION 

C KP.NO OF PARAMETERS CIJI 

DOUBLE PRECISICN T(20.20l 

DIME NSION XI 1001, Y110CI,W( 100 1, GI20I.se 1201. B 1 201,0(201. PI 10 1 

DIMENSION A(20.20l 

DIMENS IC K IPU OI 

OIMENSICN CCI20I 

ERR«0.0 

DC 1 I«1,HC 
IFIW(I > I2,2,1 

2 Wlll'1.0 

1 CCNTlHiJE 

K »0 

3 lFiK-MTI».ll,ll 

DC 5 J«1.NP 

BIJI-0.0 
DO 5 I'l.NF 



S A(l,JI«O.C 
DO 6 1-1, NO 



CALL ARG(V,XIII,O.G,PI 
OEL»(Y(ll-VI»wm 



DC t J»1,NP 
eUl^BUUDEL'DIJI 



DO 6 L-J,NP 
AIL.J1-A(L,JI«0ILI*0UI«WIII 



6 TIL. JI'AIL.JI 
P.C 7 J»l, NP 



DO 7 L»J,NP 
..V TjJ.LI'TIL.JI 



CALL INV(NP,20.T.1P.20.T) 
DO »1 J»l, KP 



CC(JI>C.O 
-D 0_M L»l, HP 



DLM»TIL,JI 
*1 CCUI«CCIJ I«eiLI»CLP 



DO U J' I, kP 
1 F1DA8SITU..II-1.CCC-30II15,1».14 



15 ERR-1.0 



|F(E RRII6.16 ,17 



17 RETURN 
16 CCNTINUE 



KCHK>0 

DC 8 J-l.NP 



IFIGU)l«14,e,914 
9H C C NTINUE 



IFI ABSICC(JI/GUII-EPSI«.8,10 
10 KCHK-1 



a GUI'GUKCCU) 
1F(KCHXU1,11^3 



11 VAR«0.0 

DO 12 I'l.NC 



CALL ARGIV.xm.D.e.PI 
___12_vaR?_V*R»>'ll l<IY( ll-VI«»2 

RGT-NO-NP 

IFIRGTIIOCI, 1001. 1002 

lOCl VAR.-O.O 

6C TO 1003 



1002 CCNTUiuE 

_V*5iS0RIil*»ZRGT) 



1003 CONTINUE 

rp 13 J-1 ,NP 



GLiM>£ABSITU,J|| 
-ij_jli-jl*>*j* SCRTICt'Mt 



RETURN 



IV C CQMPIIFH 4J£_ 



SUBROUT INE ARCIV.L.C.G.PI 



DIMENSION 0I2CI. CI20I, PIIOOI 
N»P(1I 



V'Slll 

cm-i.o 



DC 1 J-2,N 
CUI-DU-IKL 



I V'Vt-GUMCUl 
RETURN 



MEHCRY REQUIREMENTS OOOIFA BYTES 
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The Main Spectrum Analysis Program SPAN : 



\H IV C CO HPll iS.. 



LftJN., 



IM^n.^fe PA<iE_Q2Dl.„ _-'V r, COmpRFB.. 



DnijQLE PftECISICN BS(*0'>';i 

DTMEKSION NSU.?,Sni ,E(4,^ni , ASHU ,S/> I ,0» t* iSI ) ,PQLU,5n | 

DI^-fKSinN FKU,S0,2) ,tFK.5r>,ai,pr»(4,50t,DHU,S0» 
_ DIME KS ION ^>m *i 5ni.tNB» i*. i^nilCa »_i20 1 i«UAlliykL?JLiiIi^^ 

bi><Ei:sioN Dpi*,^oi,ipit(-i.i,rfi(Vo96i,ifRi*.sot .vLU(Zoi 

DfENSE'lK KllCOI.VtiaOi.ktl'lOt ,O(2n),n(20)<PUOOlf SEtZOI 
DAIfc FL,SBSE,S09F.S9SNA,SDBN»,DAV*li*VAL/ 1.0,6*0,0/ 

*i<tl,0 IVL) J)«J-1 

WRMri fe.snoii 



8001 FOftMATUm 

WRlTEt6.Z2*) 

CO ^* J"! . 

B003 FOftrtATt//, 



ZOX, 



30X,' 



TstiTOPE ANC BtANK CAT* SUKMABY 



RUN 



VPE ]ND£X= • ,14, /I 



GO TP t7T.71!,74,T*l, J 
73 REA0(S,tn2) NS(J,1.K),N£(J,2,K],EI J.KI.SMIJ 

WRITE I 6, 1023) flS(J,l,K),N^tJ,Z,K).E(J.Kt,si*ij.KI,F 



Kt J.K.I), 0F(J.K,1( 
" K,il,OF(J,>, 



132 F0RKA1t2A4, 2X. 6F10.S) 
1022 FORnATI2X.2&«,?X,2(FlC.$.2XI,2(FlCi<,<9,2Xn 

READ (5,103) PnL(j,K),CP( J,K),S[!ft( J,KJ,DSIJ.K),PaRIJ.KI,DnlJ,KI 
103 Fr.R"ATI8F)0.5) 

WRITE 16,1032) PGL {J,K) iDPtJiKI .SClPt j,K1,0S< J.K) ,PDR< JiK) , 

,1 OR ' J_lJ? ' _— 

103? FOftM7l8;"2X,Fl6.61l 



GO 1 



NSi J'i2',KY.E~(jVk)\SH(J.Kl,ASH( J.KI.DAIJ.Kt 
kT/SMi J,K )V*SHt JTk ) ,0A( J^K 1 



14 READ 1^,101) NStJ.l 
l.TAUl J, If 1,F«Y( J,K) 
MRIT£(6, 10421 NS ( Ji l.M ,NS IJ,2, 

l.TAUl J,< ).f> 1 Y{JtK) 

1042 FaRMAT(2X,2A4,7F10.5( 

104 fOflHAT12t4,?X,7F10.5l _ __ _ 

REAO (S.lb3) POL! JtKt,DP( J,Ki,SOft(J.KI,0S(J*Kl,PaR(J,K),DR<J,l 
MRITE 16,_10J2) P0llJ,K)i.OP(J,KttSOR(JtRI(OS(JjKLiPO«<JjLKI I _ 
1 OR(J,K) 

1% [F1 E(J,K1)2.2,1 _ " 

2 JPKTji=R-I 
34 COMIKUE 

WlT£i6,«57Jt 
85TJ FDRNftTllHl,///,' 

no !Ti2 j-i,4 
WRITE1&, 800 31 J 



INTEfiNAt ENERGV C*I-I6RA_TI0N._GAM| 



ITfc: 



READl!..i03l £C(J,K) 
IF(EC<J,Km76«,1764, 1963 
3 CO^I IWJE 

lTE<»,a2 1T) K, E C< J ,K1 



6217 FORHATdiX, I3,2X,F10. 

CO TC 176 3 

17t4 NCH J)»K-1 
1762 CbhTlliue 
_ _ _BEAtll%,201 ) LABt.L*ti? 
~2di FafiHtTt2A4) 

ReADl5.2C0( NFW.HOEUj 



USE.NCUT.lWD,SFH;x 



200 FaKI1iTl5T4,F10.5) 

LFM»2*_«F^M _ _ _ 

LBA5&*2*t48ASe*l 
_ LLW-2«i.NI>*I 

lDH-2*KbEL»r 
,> RlfT E(Af2U) IFWi LOEL.LSASE.NCiJT^LLWO 



,14,' 



2TTfORMATIIHi,' 

1WID?H ', 14, < SPEi 

— 24',/7//y "-— — 

____ i«T(f_<--.,6212) SFHiX 
6212 FORHATllCX, • STANDAR 

ISi^l 

Glij-C.O 

5i2J"l_,0 

bo 33 J"i, 4 
tP-IPKl jl 

lin 33 k-i, IP 

33 ITXl J,K ]-0.5* EIJ . K) 
3 REA[>(;.222IITG,NREC,< 

KDL-l _ _ _ _ __ 

NSV-N6E 
222_F0RHA1( )^,B!2,3r I0_.!i) 

HBE-HBt-il 

l_FlK(:Att3>, 31.31 

31 If (Nflta)4,<>,'5 
^ RtWlhD 4 

cotb3 " ~" 

_4 IFISBRSL U,14. IS 

1! CALL SklPia.-NSKSt"*! ' 

G O TO 3 

14 l"F(NSK»333,332"t666 
646 C*LL SK)P!0,NSK.4) 

GO TO 3 ' 

333 COHimtE 

KStt-t" '~~ 

. jH)TV-2H3'..3 34,3 37 

334 1>1NBG-A)337, 336,335 
335_Fl*t ,0 

WRUElft, 33391 
__3339 FORx.ATIZOlt. /////, 'MO FH;» CORRECTIONS!.! 
GO TO 3 
:SEI-2 



_'_!_'A!_ 



,CTY,wae,NSK,NHW,NBKS,NCAL,FTR,CaMV,TIME 



IF(ITG)6>3rfl 
R IF{LTG-ITni9S1,65atU«t 
1666 WBITE(6, 16671 
166T FORNATI^X,* TACWORD DISCREPANCY ' 

STOP__ 

65« CONTINUE 
ITN-TIM 
TMP-ITK 

FP- 10. OXTIHE-TKP 1/6,0 
T1«E-TMP*FP 
IFINCUT-NHt 9124.9124. 9125 

9124 KCUT-MCUT 
GO TO 9216 

9125 KCUT-NB 

9216 continue 

im:k-mfn*llho-kcut 

go to <711.7|21, kset 

712 ich>0.s*ecalo20..l0p0.gi 

ino'ich^zoo _ __ 

CALt BASe(LHO,r,INO.lV,BSl" " 
CAtL JM(lCH,NDELtIIf IPX! 



GO TO 717 
711 CONTINUE 

CALL DASEKWb.l.kCUT. Ifl.BS) 
IflNCALI 1_7_5* t.U54j 1755 
1755 HC-NCLI ITY) " 

10R0-H CA L«1 

PIU'tORD 
_K0L-1 . 

k«o 

DO 1756 J-1. NC 

IC-6.5*ECALtEC(ltY,Jl,LSt,GJ 

1F(IC»!WC K I40B0,40B0,]756 



CRIT''ST»RT 

I FICR IT- 10. O H 756 . 17 56, 7762 
7762" K-kti 

X ( K T . EC ( I T Y ,Ji_ 

VlKI"XMX 
, Vt-UIKl-CRIT 



«gAB15, 1031 XIKRDl.YtKROl 



1S63 CONTINUE 

L ST' LORD 

Hftif E(6,9999)" 
9999^F0ft«AT IJHl^J 



337 



6 STOP 

7 CONfiNLJE 
555 CON T INUE 



^c-' 



WRITEI6,6112) ITER, fGIJti J«lt LCRDI 
6112^ FORMAT)//, 5K,'EXTERNAtCALIBR»TI0Ni ITERATIONS • SIS.SXtTIElZ.SI 
1 ,7/ /T " ' 



9 852 f^OR'M«t(Z0](.3(F10.3.2X|,SX.F6.>,5X,r7.2l 

9BM_j:0NTINUE 

DO 3317" J-1.4 ~ ' 
IP-IPKIJI _____ 



WRITE(6,223)ITG.NREC.NBC.ITV.NSV»NSX.NRH.mK$;NCtL.FTIi;t6HV7tiNC~' 

ltiGiLC0Nt.tCON-l,tORDI 

223 FORMAT!* CONTROL ilNEl ■ 1 14. 2X.a(I 2, 1X1 * 3( 2X,Ft.4l,TI2K,FtO.»ll 

_G0 TO UTTl, 1772 1, KOt 

1772 i^RlTE(6, 17731 "' 

1TT3 F 0RKATI61K, 'ENERGY CALIBBATIOW HOT HODIFIEDM ^^_ 

1771 WRITE 16,2241 

____ _G0 TO 1718 ,7191 , KSET 

719 U~Rif~E(6,T20liCH.INX, SttC.SOD.NBE.SFLUXtFl 

— 7?P-FgRHAT..a.X..,'. FlUX CA1..I ' .■ ICH» ,. '«H*ai(.*!WC '*I4, 

15X . • -FLUX-- •,Fta.l,2X,FA.l.SX,>NBE- • .IttSK.' SFLUX> *,F10at 

2 ' FL" '.FlQ.ai 

GO TO 3 

T18 i» IT Y , 

,p^,p^,,, - 

_„__ DO 12 J -1, IP 
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SPAN List Continued 



.N IV G COMPILER 



■**l»l 



...P3T2S-LfL- 



PAfiE ODOS JH IV G COWPILtR 



GO TO IZ 
T112 CONTINUE 

CALL INI tCH|NFM,[P ilf<X) 
ItEST"! 

TFii*Bsncn-iwxi-KFxiBo n ,eoi;Taoi; 



ai! IF(FHVCIiJIJ5p07L5007i5008__ 
soot FM"v{f,jl-CbNV 
500 a C ONTmUE 



8012 KEST-: 
SOU CONTINUE 

CALL PkaStlNX.NBASE.tRteS.RATIO'.SBS.SbBI 

CRIT"SBS»l»ATIO 

GO TC It601.i6bi.41b.42nt ITV 

CU - NA CWBECTION 



Tnp-swii,j)*FMY(i,j>*eG 

81 i CONTINUE _ 

iF(SBSilOC4tlOM.i005 

1014 SBS— 0,0 

SDB — 0.0 

FK fl.J.t 1—0.0 



♦16 IF(J-1J 6,521,520 
521 SBSNA-0.P93*SeS 

SO8NA*0.'>q3*SI3a 

CO TO An 
520 iFilJ-21 *iT,9l8,*i7 
5ia SBS'SBS- SBSNA 



0F(I,J,1 J— 0.0 

pcr— 0,0 _ _ 

GO TCI 1006 ' " ----- . .... . _. 

lOOJi CONTINUE _ 

FK{i.J,i>-TbP/SBS 
CF(I.J,IHSQBTIFT«*FTR*0A[I,J1*DA1 I, JI*FKClf J.U*Fmi .J, »*S0B«SOB 

IFtrOP)3O07,30 e7,3OO fl- 

300;' PCf--0.0 "" "" '~~ '~~' ~~ "~ 

GO TO_100^_ . ■ 

30011 PCt-lbO.«Bfi/TOP ' " 

lOOfi CONTIMJE 



S08'5OHT(SOB*S0B»SDeNA*SOaNA) 

sesN*-f>.o 

SDSNA-O.O 
GO TO 417 
HG - SE COHBECTION 
42 1 1F(J~ 11 <Oq,*OT,40q 



407 Sa$E-a.185*S6S 

SDBE-0.3a5«SOB 

GO TO 417 
409 If(J-3l 417,4^8,417 _ 
40B SflS-SBS-SBSE 

SBo-SQftrisap»SR9tsqBE*spae' 



SBSE-O.O 

SOBE-O.'l 
417 NBE-1 

GO TCr415,415.415.412.412}, NBG 
415 NB-^NBRd.JI 

_G0 TO ( 4091 , 40921 , Nfl 



wBitEi6,!sn"2) "CRiT ■■ " .^~- _— _ 

.92i>_f OftNATdXtl 3t2X, 14, IX tLti-SX f 2**. lX.Fe.Z, 1X, F9. I. IK,F6. I , 

l5X,El2.5.2)t,El2.5,5X."F6.2I 
GO TC 12 _^_______ ■ 



4091 



CONT INUE 
ARG-CTIME-TVi 

ifiabsiaSci-. 

1001 FCOB-1.0 

GO TO 1003 
10Q2 FC0«-EXP<ARG1 



8014, CONTINUE '■ 

_ . WBITEr6,4Q04tJ.ICH.im.WSII .I.J> ,HS 11.2.JI,EH,JI,SBS,SDB, LFM 
MRITE(6,57I2» CBlt - - ■ 

N8Rn,Jl-2 

i; CONTINUE ' ~^ — _ - — — _ 

_ WBITE 16.29481 



1003 GO TO <4a01, 40021, KEST 

4001 CONTINUE _ 

TF-FKi I,J,NBEl«fL*FCDB 

TD-Of (1,J,NBE(*FL*FC_PR 

TOT-SBSilF 

T0DS-SBS*S8S*TD*T0*Tf»TF»SQB«SD a 



URITE(6,4941) IL VLIJl , J- 1.231 

*'941 FORHATI/, 12X. 2 3112 , 3 )(l,/ t _ 

00 2947 J-I. IP 
IFHIPIJH75I4.75I4. 2947 



-l/l 



GO TO (413,414,815). N6G 

413 6G-FTft*P0L(l,J) __ 

ANET-TOr-BG'" "" - - - 

ONET-S0RT(TODS*FTR*FTP*CP( b •* »*i"!l I f JJJ 
GO TO 420" '" ' ' ■"■ 

414 BG-FTR*S0R 1 I.J1 



7514 CONTINUE 

IST-ITR(I,JX-11 _ 

INb-ISf*22 '—^ 

HRITE(6.J949( J, 1ST, ( IR(L 1 ,L« 1ST, INDI 

2947 CONTINUE ' " " 

a 949 F0ftMAT(lX.[2,lX,I4,2X,23I51 



GO TO 3 
_224 F0RH AT1//1 
ENO 



■ t MEMORY BEQUIBEMENTS 013768 BYTES " 



(l.JI) 



Cn TO 4 20 
815 BG'-FTR^PPfHI.JI 

ANET-T0T-6G 

bNEt-SgRT(T6DS»FTP»FTS»DHri,"j)»DR(i 
42 _V AL -AN ET /CONV 



-SJ^BRQUnNE-XmiHX. IB. XMX.NPf. I 



DVAt-ONET/CONV 
; AL - HG CORRECTION 

ifliTy-it 605.666,605" 

606 IF(J-3) 615,607,615 

607 AVAL-a.2*VAL 
CAVAL-Q. 2»0VA L 



GO TO 615 
605 IFlITY-2) 615,608,615 
60(! IF{J-6» 615, 669,615"' 
609 VAL-VAL-AVAL 

OVALiSOHT(bVAL»DVAi*0AV* 
CAVAL'O. O _ _ 

AVAl-b.O 
615 IFITOTl 2001,2001,2002 

zooi Pct--o.b " 

GO TC 2003 
2602 CONTINUE ' 

PCT-100,»BG/TGT 



ti:p-ibu3}-ir(ji) 
. B ("Jl- liLair^IAl J PJU .t.I Rl Jil_ 



lr- tBQT Il,2,2 



2 NI'C-l 
J<t'Jt-0.0_ 

3 fifruRN 



"HfcIlCM.. R£flUll^£MEMI^ CC029ft (} YT£5 



NS{I,l,JI,NS(I.2,J),Et ItJI.SBS.SDB.VAL 



2003 CONTINUE 

WfirTE(6,500l J.ICH 
I ,DVAL,PCT,FCOB 
500 F0HHAT(1X,I3,2X,K,1X.I4,_3X,2A4,IX,F8.Z.1X,F9.1,1X,F6.I, 
1 5X,F11.4,IX,F11.4.2X,F6.ir6X,F6.4l " '" 

WR1TE( 6, 5 7 121 CHI T _-___^ 

5712 F0fit»AT{lH*,115X,Fl0.21 "~ ~~ 

JFtNRW11155, 1157, 1157 
1155 CONTINUE' - --- — _ __ _ 



Si)BRi3U.UN£. _U!Ll£KiflS£j.l8Jflxj_ 



DINENSICh IB(4096) 

"!;-iPK -Nsc 

01PK4SSC 



_ Ji t6_f 0«"*I' 2*4, 2X,3F10. 2, 2flO_^4^X, 2 A4, IX, *4) 

1157 CQNtiNUE " ' "" ~~ 

GO TO 12 
4002 W«rTE16,4004U,ICK,IMX,NS(r,l,JUNS(lVZ,JUEnVjI,SBS,SDB,LFw" 
•4004 F0RMATIU,I3,2X.I4,1X,I4,3X.2A4,1X,F8.2,1X,F9.1,JX.F6.1, 



1' PEAX NOT WittilN 

"?. 1 1 UJL' S 712) CRI T 
CO TC 12" 



,14,' CHANNEL SCAN "iNTERVAt 



- fe£H£M_REQiiEREMENTS Qor 2QE BWES 



4092 CONTINUE 

WBnE(6,"40?lVjVlCH,IMX,NSn;i;jUNSl'i.2.ji,y(T.7r,TBS'.V0B~ 
4021 f0'"'ATfIX.,13,2X.I4,lX,I4,3X,2A4,lX,F8.?,lX,F9,Ui)(.FAaJ 

!• CALC. DELETED - NC PEAK IN STD. ', - t.iA_-._-_tAJ.__ 

H«1TE1 6,5T 12I CRIT _„__^_ 

GO TO 12 ' '■ ■"" ' 

412 CONTINUE 

GO TO (a0i3,8014l'.' KESt ~ "" 

8013 CONT IMJE 

NBR(I,J)"1 " ~"" ' — 



iw c cnttPti pp 



-_ fU«.IlQN_££»lUjMAl^ei 

DlMENSICh Gf20I 

X-liC __ 

EC*L-"G(1I 
op 1 -J-Z. hCAL 

X-xSe ' 
\ EC AL-ECAL*C(JI» x 

RETURN ~~ ~ 
ENO 



